Attractor scenarios and superluminal signals in /c-essence cosmology 
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Cosmological scenarios with fc-essence are invoked in order to explain the observed late-time 
acceleration of the universe. These scenarios avoid the need for fine-tuned initial conditions (the 
"coincidence problem") because of the attractor-like dynamics of the fc-essence field (j>. It was recently 
shown that all fc-essence scenarios with Lagrangians p = L(X)(j>~ 2 , where X = ^4>, IJ ,(f>' >i , necessarily 
involve an epoch where perturbations of (j> propagate faster than light (the "no-go theorem"). We 
carry out a comprehensive study of attractor-like cosmological solutions ("trackers") involving a 
fc-essence scalar field <j> and another matter component. The result of this study is a complete 
classification of fc-essence Lagrangians that admit asymptotically stable tracking solutions, among 
all Lagrangians of the form p = K(4>)L(X). Using this classification, we select the class of models 
that describe the late-time acceleration and avoid the coincidence problem through the tracking 
mechanism. An analogous "no-go theorem" still holds for this class of models, indicating the existence 
of a superluminal epoch. In the context of fc-essence cosmology, the superluminal epoch does not 
lead to causality violations. We discuss the implications of superluminal signal propagation for 
possible causality violations in Lorentz-invariant field theories. 



I. INTRODUCTION AND OVERVIEW OF 
RESULTS 

Cosmological scenarios involving a scalar field known 
as fc-essence [l|, 0, H[ are intended to explain the late-time 
acceleration of the universe (see Ref. |4| for a recent re- 
view of dynamical models of dark energy) . An important 
motivation behind the fc-essence scenarios is to avoid the 
fine-tuning of the initial conditions for the scalar field 
(the "coincidence problem"). 

The effective Lagrangian p(X, (f>) describing the dy- 
namics of the scalar field <fi consists of a noncanonical 
kinetic term, 



p{X,4>)=K{ct>)L{X), X 



1 



:0 L 



(1) 



where K((f>) and L(X) are functions determined by the 
underlying fundamental theory. One considers the evo- 
lution of the field 4> coupled to gravity in a standard ho- 
mogeneous cosmology in the presence of matter. With a 
suitable choice of the Lagrangian, the evolution of <p dur- 
ing radiation domination quickly drives the system into 
a region in phase space where the fc-essence field <j) has 
a nearly constant equation of state with = i, mim- 
icking radiation. Thus the energy density of fc-essence 
approaches a constant fraction of the energy density e m 
of the radiation. This behavior of fc-essence (w$ — > const 
and e^/etot — * const, where £tot = £<p + £ ro ) is called 
tracking, and the solution with const is called a 

tracker solution. 

The parameters of the Lagrangian can be adjusted such 
that the energy density in fc-essence during the radiation 
era is small (e m ~ £tot), so that the standard cosmo- 
logical evolution is not significantly altered. After the 
onset of dust domination [w m = 0), the energy density 
in fc-essence quickly becomes negligible and the evolu- 
tion leaves the radiation tracker. A tracking solution 



with = does not exist (due to a particular choice of 
the Lagrangian) , and instead the fc-essence is driven to a 
tracking regime with w^, ~ const < 0. Since < w m , 
the fc-essence will eventually dominate the energy den- 
sity of the dust component. The precise value of in 
that regime can be paramctrically adjusted to fit the cur- 
rently observed data; in particular, values m^, ~ -1 can 
be achieved. 1 

In our terminology, a "tracking solution" is a solution 
for which approaches a fixed value, whether or not 
this value is equal to the equation of state parameter w m 
of the dominant matter component. It is essential that 
the tracker solutions are stable attractors for nearby solu- 
tions. Because of this property, the field <fi is driven into 
the tracker regime in the phase space with fixed values 
of Wrj> and the ratio e^/etot, for a wide range of initial 
conditions for <fi. To construct a viable fc-essence model, 
it is important to choose a Lagrangian p{X, <f>) for which 
stable tracker solutions exist within the radiation- and 
dust-dominated cosmological eras. 

Previous works concerning the dynamics of fc-essence 
either assumed a specific form of the Lagrangian, for in- 
stance Q 



p{X,<t>) 



L(X) 



(2) 



or imposed ad hoc restrictions on the Lagrangian with 
the purpose of deriving exact solutions (e.g. P). In 
particular, it was assumed that = const is an ex- 
act solution of the equations of motion. However, the 



1 We note that the "phantom" values w 
in this single-field model; see e.g. 0, 



< — 1 cannot be reached 
"Phantom" models such 
as that of Ref. Q cannot describe the tracking behavior of fc- 
essence since in these models iuj, < — 1 at all times. 
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physically necessary requirement is weaker: namely, one 
merely needs that should approach a constant value 
asymptotically at late times. The existence of an exact 
solution uirj, — const is not necessary. With this weaker 
requirement, a much wider range of Lagrangians enters 
the consideration. 

In the present paper, we restrict our attention to La- 
grangians of the "factorized" form (|TJ) but do not impose 
any further a priori restrictions on the Lagrangians; nei- 
ther do we require the existence of analytic exact solu- 
tions, or of solutions with — const. It is only as- 
sumed that the cosmological scenario is realized with 
4> > and that <f) reaches arbitrarily large values. Our 
results can be viewed as a comprehensive extension of 
previous studies of attractor behavior in fc-essence cos- 
mology (e.g. d, [^, [HI). We determine the class of La- 
grangians p(X, (j>) that admit stable tracking regimes in 
which u>0 — * const, for a given value of w m . The possible 
asymptotic values of w$ and e^/etot are derived in each 
case. 

The form ([T]) is sufficiently general to reproduce an ob- 
servationally measured cosmological history [Tlj and cov- 
ers many interesting cases, such as fc-essence with purely 
kinetic term [l2j or the "kinetic quintessence" [l[ . Factor- 
ized Lagrangians have been the main focus of attention 
in the study of fc-essence (see e.g. H d). More 
generally, Lagrangians of the form 

p(X, <j>) = [KxiflX** + K 2 {^)X^r , (3) 

where 711,712,713 are constants, can be reduced to the 
Lagrangian |T]) by a suitable redefinition of the field </>. 
Our analysis will also apply to Lagrangians that have 
the asymptotic form p m K((f>)L(X) for 4> — > 00 and 
for which only the large- <p regime is cosmologically rele- 
vant. Nonfactorizable Lagrangians, such as those studied 
in Refs. [1, [H, [l7], d, [Hi, require a separate considera- 
tion which we do not attempt here. 

Recently, it was shown that the scenarios of fc-essence 
cosmology with Lagrangians of the form @ neces- 
sarily include an epoch when perturbations in the fc- 
essence field propagate faster than light (the "no-go the- 
orem" [20]). It is well known that super luminal propaga- 
tion of perturbations opens the possibility of causality vi- 
olations, although causality is actually preserved in many 
cases. This issue has been a subject of some debate, see 
e.g. the discussion in Refs. [HI H2, HI II H, H M, HI • 
One of the motivations for the present work is to de- 
termine whether the "no-go theorem," derived for a re- 
stricted class of fc-essence Lagrangians, still holds in sce- 
narios with more general Lagrangians. 

To answer this question, we performed an exhaustive 
analysis of all the possibilities for the existence of stable 
tracking solutions in ghost-free fc-essence theories with 
positive energy density (the complete list of physical re- 
strictions is given in Sec. [TTJ) . We considered the cos- 
mological evolution of a scalar fc-essence field <fi coupled 
through gravity to a matter component having a fixed 
equation-of state parameter w m . In this context, we enu- 



merated all Lagrangians of the form (TTJ) that admit at- 
tractor solutions with w,p — ► const and e^/stot —* const 
at late times (Sec. IV Aj) . Since our task is to determine 
the entire class of theories admitting a certain asymptotic 
behavior, numerical calculations could not be used. The 
analytic method used for the asymptotic analysis of the 
dynamical evolution is outlined at the beginning of Ap- 
pendix \X\ where all the calculations are presented in de- 
tail. This method is similar to that developed in Ref. (29[ 
for the analysis of attractors in models of fc-inflation. 

Armed with the complete enumeration of stable track- 
ers, we then selected the Lagrangians capable of provid- 
ing a subdominant tracker solution during the radiation 
era and an asymptotically dominant tracker solution dur- 
ing the dust era. We show that the only appropriate class 
of Lagrangians consists of functions p{X, 4>) of the form 

p(X, 0) = 1 + f 2 ° W L(X), lim K (4>) = 0. (4) 

Since the dynamical evolution drives cj> towards very large 
values, these Lagrangians are practically indistinguish- 
able from the Lagrangians of the form ((2]). Then one can 
prove, similarly to Ref. [20], that the cosmological evo- 
lution necessarily includes an epoch where perturbations 
of the fc-essence field <j) propagate with a superluminal 
speed. Thus, we prove the "no-go theorem" starting from 
a much wider initial class of fc-essence Lagrangians. 

In Sec. IIIII we discuss the implications of superlumi- 
nal signal propagation for causality. The cosmological 
scenario of fc-essence does not exhibit any causality vio- 
lations at the classical level, despite the presence of su- 
perluminal signals. Preservation of causality in a general 
configuration of fc-essence field can be viewed as a po- 
tential problem, on the same footing as the chronology 
protection problem in General Relativity [30] . 

II. PHYSICAL RESTRICTIONS ON 
LAGRANGIANS AND SOLUTIONS 

In this section we consider some physically necessary 
restrictions on the possible Lagrangians p(X, <p) and so- 
lutions <j>(t). 

The main physical context for fc-essence scenarios is 
the evolution of the fc-essence field on the background of a 
matter component with a fixed equation of state parame- 
ter w rn . The energy density of fc-essence is not necessarily 
dominant during this evolution. Since fc-essence scenar- 
ios are proposed as an explanation of the dark energy, 
we do not consider the case w m = — 1 (during primordial 
inflation, one must also have w m > — 1 due to the ne- 
cessity of the graceful exit). However, we leave open the 
possibility w m < — 1. 

An important requirement for a field theory is stabil- 
ity. A theory for a field <j> is stable and ghost-free if the 
energy density is positive, the speed of sound c s is 
real (not imaginary), i.e. c 2 s > 0, and the Lagrangian 
for linear perturbations has a hyperbolic signature and a 
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positive sign at the kinetic term. The speed of sound for 
perturbations on a given background is given by 31 1 



III. SUPERLUMINAL SIGNALS AND 
CAUSALITY 



P.v 
Vp.vi 



(5) 



To obtain the leading terms of the Lagrangian for the 
perturbations, one writes a perturbed solution as <f> = 
(f>o(t) + x(^: x ) an d expands the Lagrangian p(X, <fi) to 
second order in \] the Lagrangian p(X, <fi) is assumed to 
be an analytic function of X at X = 0. The relevant 
terms are those quadratic in the derivatives of x-i 



p{X,(j>) 



1 



P.XTpC^X 1 " + \^P,XXX^X,vK4>l + ■■ 

^•• ; "'\,.\, • ... 



(0) 



It follows that linear perturbations x propagate in the 
effective metric 



G^ v = p,x9^ v + P,xx<t 



(7) 



The no-ghost requirement is that the metric G^ v should 
have the same signature as g^ v . Regardless of whether 
the 4-gradient d>£ is spacelike or timelike, 2 the resulting 
conditions are [32j 

P,x = ~P,v > 0, p.x + 2Xp.xx = P,w > 0. (8) 

v 

In the cosmological context, the field is a function of 
time t only; in standard fc-essence scenarios that we are 
presently considering, <j)(t) grows monotonically with t. 
Hence, <^>g is timelike and the velocity v = (j) is positive, 



v 



do 



(9) 



We conclude that a physically reasonable cosmological 
solution should satisfy (for v > 0) the conditions 



vp, v - P > 0, p :V > 0, p vv > 0. 



(10) 



It follows that p(v, (f>) v= o < and that p(v, 4>) is a con- 
vex, monotonically growing function of v at fixed (j) (&t 
least for values of <fi and v relevant in a cosmological sce- 
nario). For factorized Lagrangians p(v, <j)) = K((f>)Q(v), 
we find that Q(v) must be a convex, monotonically grow- 
ing function of v with Q(0) < 0, and also Q'(v) > and 
Q"(v) > for all values of v > that are relevant in a 
given cosmological scenario. 

Finally, we assume that K{<p) has monotonic behavior 
at d> — > oo . 



One of the results of this work is a conclusion that ev- 
ery fc-essence scenario based on attractor behavior and a 
Lagrangian of the form (fT]) will include an epoch where 
the perturbations of the fc-essence field propagate super- 
luminally. It is therefore pertinent to discuss the possibil- 
ity of causality violations in the presence of superluminal 
signals. 

We first consider small perturbations (f>o + 8<f> of an ar- 
bitrary background solution 4>q(x) in a Lorentz-invariant, 
nonlinear field theory. To first order, the evolution of S(j> 
is described by a linear equation of the form 

G^cfo^^M + £"[0o]V M <ty + C[0 O ](50 = 0, (11) 

where the coefficients G^ , B^ L , and C are determined by 
the Lagrangian and depend on the background solution 
(j>o . Unless Eq. (|11[) is hyperbolic (the matrix C" having 

signature H or equivalent), the theory will trivially 

violate causality: an initial-value (Cauchy) problem will 
be ill-posed in any reference frame, and the evolution of 
perturbations will be physically unpredictable. There- 
fore, it is necessary to require that G^ v have a hyper- 
bolic signature. Background solutions <po that lead to a 
parabolic or an elliptic signature of G^ v even in a small 
spacetime domain must be avoided as pathological. The 
cosmological solution (t>o{t) used in fc-essence scenarios 
will be well-behaved if the conditions (jTUJ) hold. Below 

we assume that G^ v has signature H . 

Within a sufficiently small spacetime domain, we may 
regard G^ v , , and C as constants. Then it is straight- 
forward to derive the dispersion relation 



G^ h^ky 



iB^k,, 



C = 



(12) 



for plane wave perturbations 8<j){x) cx exp [ifc^x^]. In or- 
der to send information ("signals" or "sounds") by means 
of a perturbation S(j)(x), one needs to create a wave front, 
that is, a perturbation with an extremely short wave- 
length and a high frequency. Thus, wave fronts prop- 
agate along wave vectors fc^ determined by the leading 
term in Eq. (TT2]), 



G^k^h = 0. 



(13) 



Any wave packet consisting of a superposition of plane 
waves will propagate behind the wave front. Therefore, 
a 4-vector u M of signal velocity must lie within the sound 
cone, 



G^UftUv > 0. 



(14) 



Since the "sound metric" G^ v is determined by the lo- 
_ cal behavior of the background solution 4>q{x), the sound 

cone may have an arbitrary relationship with the light- 
is null, the metric will have the correct signature if con e g^Uy.Uu = determined by the spacetime metric 
> 0. g 1 *" . Thus, in some theories the sound signal worldlincs 



If , 

p x > and p 
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may be timelike, null, or even spacelike depending on the 
spatial direction of their propagation. 

The speed of sound waves is therefore direction- 
dependent. The background tensor G^ v determines (a 
class of) preferred reference frames where G^ is diago- 
nal. Propagation of sound is most conveniently described 
in terms of sound speeds in different directions in a pre- 
ferred frame. In this sense, one may say that a dynamical 
Lorentz violation takes place for sound waves, although 
the full theory (that includes the tensor G^ v as a dynam- 
ical quantity) of course remains Lorentz- invariant. 

In the context of fc-essence cosmology, the sound metric 
Q^v i s given by Eq. ([7]) . Preferred frames are those where 
the t axis coincides with the cosmological time. So in the 
preferred frames <po — (f>o(t) is a function of time only, 
and the dispersion relation is 

^ 2 =^|k| 2 , (15) 

where k is the 3-dimensional wave vector and c s is 
the (direction-independent) speed of sound defined by 
Eq. ([5]). In this paper we show that the considerations of 
the "no-go theorem" [2Gj hold for those Lagrangians of the 
form ([I} that admit scenarios of tracking fc-essence. By 
virtue of this theorem, there exists an epoch with c 2 > 1. 
During this epoch (which may be quite short 20}), it is 
possible to send signals along spacelike worldlincs. 

If spacelike sound signals propagated in arbitrary 
spacetime directions, one could easily create closed world- 
lines made of signals, called "closed signal curves" (CSCs) 
in Ref . [28| . This would open a Pandora's box of classi- 
cal time travel paradoxes, also violating the unitarity of 
quantum theory (see e.g. [H, [H], HH, HH ) . However, the 
allowed sound signal directions are only those within the 
sound cone (1141) . This limitation precludes the possibil- 
ity of constructing CSCs within a small domain where 
G^ v ~ const. This can be shown as follows. Diago- 
nalizing the tensor G^ v within that domain, one finds a 
preferred reference frame {t, x, y, z} where sound signals 
(whether spacelike, null, or timelike) always propagate 
in the positive direction along the t axis. Signals sent 
by conventional means also propagate in the positive t 
direction. Since the local coordinates {t, x, y, z} are valid 
within the entire domain where « const, no CSCs 
are possible within that domain. 

It is straightforward to see that no causality viola- 
tions through CSCs can occur in fc-essence cosmology. 
Since 0o > at all times, the 4- vector V^o is every- 
where timelike and selects a global preferred reference 
frame. (Even if 0o — momentarily, the preferred frame 
is still selected by continuity.) In this reference frame, 
the sound waves propagate in the direction of increasing 
coordinate t. Hence, there exists a global foliation of the 
entire spacetime by spacelike hypersurfaces of equal t. 
Any sound signals (whether spacelike, null, or timelike), 
as well as any signals sent by conventional means, will 
traverse these hypersurfaces in the direction of increas- 
ing t. It follows that CSCs cannot occur, either locally 
or globally. 



Similar conclusions were reached in models of inflation 
having c 2 > 1 as well as in situations involving a 
fc-essence field on a black hole background [23], ■ By 
itself, a superluminal speed of sound does not automati- 
cally lead to CSCs or causality violations. 

In certain field theories, one can construct backgrounds 
4>o(x) where CSCs are possible; a notable example is 
given in Ref. (22j. However, such backgrounds are ar- 
tificial in the sense that they require an ad hoc configu- 
ration of the field <fio(x). It remains to be seen whether 
such causality-violating backgrounds can occur as a re- 
sult of the dynamical evolution of the field <fio(x) in a 
cosmological context. 

The problem of causality violation by CSCs is similar 
to the problem of closed timelike curves (CTCs) occur- 
ring in General Relativity [3fJ ■ It is difficult to find a met- 
ric g^u that is initially well-behaved but admits CTCs as 
a result of dynamical evolution (one such example is given 
in Ref. [38(). Hawking's "chronology protection conjec- 
ture" states that such spacetimes containing CTCs will 
be always unstable due to quantum effects; but it remains 
an open conjecture [30| . Similar considerations apply to 
CSCs occurring in nonlinear field theories. It is possible 
that CSCs will always lead to quantum instabilities due 
to a similar "chronology protection" mechanism. Further 
work is needed to resolve this intriguing question. 

IV. EQUATIONS OF MOTION 

We begin by writing the well-known evolution equa- 
tions for fc-essence cosmology in a convenient set of vari- 
ables. The equations in this section will be used at vari- 
ous points in the following analysis. 

We consider a spatially flat FRW universe with the 
metric 

g^dx^dx" = dt 2 - a 2 {t) [dx 2 + dy 2 + dz 2 ] , (16) 

where a(t) is the scale factor. In the epoch of interest, the 
universe contains the dynamical fc-essence field <p{t) and 
a matter component with energy density s m and pressure 
p rn . The matter component can be approximately treated 
as nondynamical in the sense that its equation of state is 
fixed, 

w m = = const. (17) 

The energy-momentum tensor of the field 4> is that of a 
perfect fluid with pressure p{X, </>) and energy density 

e = 2Xp, x - p. (18) 

Here and below we denote partial derivatives by a 
comma, so p.x = dp/dX. We introduce the velocity 
v = <f> as shown by Eq. ©. Note that the Lagrangian 
p(X, <t>) is an analytic function of X and thus an analytic 
function of v 2 . 
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The equation of state parameter for fc-essence, w^, is 
defined by 



W4, 



p{X,<t>) 

£(/> 



P 



Vp,v - P 



(19) 



A factorizable Lagrangian (TTJ) is expressed as a function 
of v and (ft as follows, 

p(X, <f) = K(cft)Q{v), Q(v) = L{X). (20) 

For a Lagrangian of this form, is a function of v only, 

Q 



w<f,(v) 



(21) 



vQ'-Q' 

since the energy density factorizes, 

e$ = K(<f>)i 4> (v), e^{v)=vQ'{v)-Q{v). (22) 

We assume that the functions K and Q in Eq. (|20|) are 
chosen such that -ftT(</>) > 0. 

The cosmological evolution is described by the equa- 
tions of motion for 4>{t), £ m (t), and a(t), 



- = H = K^/e^ + e,„, 



8ttG 



(23) 



— (p, v (v,4>)) 



>P,w + <j>P,<jw = ~3Hp iV + p.tp, (24) 

-3H (e m + p m ) = -3H (1 + w m ) s m . 

(25) 



The equation of motion for the field (ft can be also rewrit- 
ten as a conservation law, 

e = -3H (e + p{X, d>)) = -3H (1 + W(j> ) e . (26) 

The total energy density e t ot = £<£+£m satisfies the equa- 
tion 



£tot — — 3i?£t 



(1 + w rn ) + (w<i, ~ w m ) 

Etot 



(27) 



Since the equations of motion (|2"3")) - (r2"5"|) do not depend 
explicitly on time, and since (ft{t) is monotonic in t, we 
may use the value of (ft as the time variable instead of t. 
Then we obtain a closed system of two first-order equa- 
tions for v{4>) and e m {(ft), 

dv{4>) _ vp, vt f, - p t <f, + 3np, vy /e m + vp :V - p 



vp 



de. 



3k (1 + w m ) e r 



(28) 

y/e m + vp, v - p. (29) 



We will make extensive use of the auxiliary quantity R 
defined by 



I? 



(30) 



Since energy densities e§ and e m are always positive, the 
ratio R always remains between and 1. The equa- 
tion of motion for R(4>) is straightforwardly derived from 
Eqs. (|25p - (|2"6"|) and can be written as 

S = -— R{l-R){w m -w <t> {v, ( j>)). (31) 
do v 



We may reformulate the equations of motion (|28|) (129)) as 
a closed system of equations involving only the variables 
v((ft) and R(4>). Since 



1 - R 



VP,v ~ P 
l~R ' 



we obtain 
dv 1 



vp, vv 



vp,v4> - P,4> + 3np,v 



dR 



3 k 



— = R\/l - R^Vp^y - P 



VP.v - P 

1-R 
P 



VP,v ~ P 



(32) 



(33) 



(34) 



For Lagrangians of the form (|20| . these equations are 
rewritten as 



dv 2 
12 = - c s(v) 



1 + w^v) 



1 - R 



(35) 



dR 



—RVi^R^ K((t>)£<t,(v) (w m - wt(v)) . (36) 



Here e^v), c 2 s (v), and w^(v) are understood as fixed 
functions of i>, 



e^v) ee vQ' - 



(37) 

determined by the given Lagrangian p(v, (ft) — Q(v)K((ft). 
These functions satisfy the following equations, 



d _ 1 + w (w) „ 

dv vcf{v) 

d 1 + wJv) \ wAv) 

— w^(v) = * 1 f-— 

dv v c i( v ) 



(38) 
(39) 



V. VIABLE LAGRANGIANS FOR TRACKING 
SOLUTIONS 

The detailed analysis of asymptotically stable solutions 
is given in Appendix [Al Each asymptotically stable so- 
lution is characterized by the asymptotic values of v = (ft 
and R — e m /etot, considered as functions of (ft: 

vo = lim v(cj>), i? ee lim (40) 

<j>-KX 0-+oo £tot{(ft) 

As a summary of the results, we list all of the possibil- 
ities, together with the requirements on the Lagrangian 
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p = K((f>)Q(v) and the allowed values of vq, Rq, w m , 
and w^Vq). [Note that the function K{<p) can be al- 
ways multiplied by a constant, to be absorbed in Q(v).] 
The requirements listed are necessary and sufficient con- 
ditions for the asymptotic stability of tracker solutions. 
The applicability of these tracker scenarios to fc-essence 
cosmology is analyzed in subsections IV Bl IV CI and IV Dl 



Tracker solutions 



Case 1. The function K(<fi) is of the form 

= Kq{4)) = Q (41) 

0^ <f>— >oo 

The value vq is determined from w^vq) — w m , and then 
Ro is given by 



9k 2 Q'(v ) 2 
4£«/,(wo) 



Ro = 1 



(42) 



This value of Rq must satisfy < Rq < 1 (the possibility 
Ro = is equivalent to case 2). The conditions 

vq ^ 0, c 2 s (v ) ^ 0, \w m \ < 1, w m < c 2 (v ), E^Vq) ^ 

(43) 

must hold. 

Case 2. The function K((f>) is of the form 

= lim Ko{4>) = (44) 

0^ 0^OO 



The value wq is determined from 



V 1 + W0(w o ) 



and must satisfy vq ^ 0. The following conditions must 
hold, 

u^(uo) < w mi \w^,{v )\ < 1, c 2 (v ) ^ 0. (46) 

The tracker solution has Rq = (fc-essence dominates at 
late times). 

Case 3. The function K(<fi) is of the form 

K(*)-™ Bm^-D, ,47) 

i.e. the function Ko either tends to a constant, or grows 
or decays slower than any power of at — > oo. This 
condition determines the value of a. This value of a must 
satisfy 



2 < a < 1 



1 + w n 

The interval for a is nonempty if 
\w m \ < 1. 



(48) 



(49) 



(50) 



The value of a determines vq by 

01 — Z 

1 + w m 

The resulting value of vq must satisfy the conditions 

v ^ 0, c^(u ) > Wj,(v ) > w m ,c 2 (v ) ^ 0, i^(v ) ^ 0. 

(51) 

The tracker solution has Ro = 1 (fc-essence is negligible). 
Case 4- The function if(^) is of the form 



K{<j>) 



(52) 



where the function Kq must satisfy 



lim K Q (4>) = 0, lim ln f° ( ^ = 0, (53) 

i.e. K ((j)) decays slower than any power of at — > oo. 
The value of vq is determined from the conditions 

Mty(«o) = Wm-, \w rn \ < 1. (54) 

The following conditions must then hold, 

v ^ 0, c^(u ) > w m , e^vq) ^ 0, Cg(u ) 7^ 0. (55) 

The tracker solution has Ro = 1 (fc-essence is negligible). 

Case 5. The function K(<f>) decays slower than 0~ Q 
(or grows), where 



1 + w„ 



-1 < w m < 0. 



(56) 



(45) More precisely, 



= lim K o(4>) = oo. (57) 



The value of vq is determined as a root of Q(vo) = and 
Q'(vo) = 0, i.e. we must have a Taylor expansion near 
v = vo of the form 

Q(v) = —{v-v ) n , n>2, Q o >0. (58) 

Then the tracker solution has Ro — 1 (fc-essence is negli- 
gible) and iv<f,(vo ) = 0. The value of w must be above vq 
at all times (or else c 2 < 0). 

Case 6. The function Q(v) has an expansion at v = 
of the form 

Q(«) = Qi«" + o(v n ), Qi>0, n>2. (59) 

This determines the value of n. The function if (0) de- 
cays slower than 0~ Q (or grows), where 



2n 



(n- 1) (1 + w m ) 



(60) 
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More precisely, 



K(4>) = i™ k (4>) = oo. (61) 



The condition 



n- 3 1 
r < w m < r 

71—1 71 — 1 



(62) 



must hold. Then the tracker solution has Rq = 1 (fc- 
essence is negligible), vo — 0, w^{vq) = -^hj, and c 2 (vq) = 
l 

n-l ' 

Case 7. The function K{<f>) has the form 



K{cj>) 



Ko(<t>) 



(63) 



where the function Kq{4>) is such that 
1 



lim K {4>) > 



or lim Kq{4>) — oo. (64) 



4>— >oc * ' 9ft Qi 0— *oo 

The function Q(w) has an expansion at v = of the form 

Q{v)=Q 1 v 2 + o(v 2 ) 7 Qi >0. (65) 

We must have w m > 1. The tracker solution has i?o = 
(fc-essence dominates), t>o = 0, and w<j>(vo) = c 2 (vq) = 1. 

Case 8. The function K(cj)) decays slower than <\>~ 2 
or grows, 

= lim K °(<t>) = oo- (66) 

The value of i>o is determined from C/'(t>o) = 0, Q{vq) < 0. 
More precisely, we have an expansion near v = vq, 

Q{v) = Q Q + Q 2 {v-v Q y\ Q o <0, n>2. (67) 

We must have vq ^ and w m > — 1 . The tracker solution 
has Rq = (fc-essence dominates) and w</,(fo) = —1- The 
value of v must be above vo at all times (or else c 2 s < 0). 

Case 9. The function K(cf>) decays slower than <\>~ 2 
or grows, 

K(4>) = lim K (4>) = oo. (68) 

The value of vq is determined from Q'(i>o) = 0, Q(vq) = 0. 
More precisely, we have an expansion near v = vq, 

Q(v) = Q 1 (v-v ) n , n>2. (69) 

We must have w rn > and vo ^ 0. The tracker solution 
has i?o = (fc-essence dominates) and w^,{vq) = 0. The 
value of v must be above vo at all times (or else c 2 < 0). 

Case 10. The function K{<f) decays slower than <f)~ 2 
or grows, 

K Q {<t>) 



K{4>) = 



lim K $((/)) = oo. 

0— >OC 



(70) 



The function Q{v) must have an expansion near v = of 
the form 



Q(«) = -Qo + Qi^", Qo>0, ri > 2. 



(71) 



We must have w m > — 1. The tracker solution has Rq = 
(fc-essence dominates), t>o = 0, and w^vq) = — 1. 



Case JJ. The function X(</>) decays slower than (f> 2 
or grows, 

K{4>) = lim K (4>) = oo. (72) 

The function Q(i>) must have an expansion near v = of 
the form 

Q(v) = Qiv n + o(v n ), Qi>0, n>2. (73) 

This determines the value of n. The condition 

1 



n- 1 



(74) 



must hold. The tracker solution has i?o = (fc-essence 
dominates), vq — 0, and w^,(vo) — c 2 — zrzj. 

Case 12. The function Q(i>) must have an expansion 
near v — of the form 

Q(v) = Q 1 v n + Q 2 v n+ P, Qi > 0, n > 2, p > 0. (75) 

This determines the values of n and p. The function K(<p) 
must be of the form 



where Ko((j)) must satisfy 



(76) 



lim K (<j>) = oo, / -ZK "- 2 (0) = oo. (77) 

<£->°o J (p 

(The function Ko(cj)) grows slower than (\n(p)^ n 2 ^ p .) 
We must have w m = . The tracker solution has Rq = 

(fc-essence dominates), vq — 0, and ^^(wo) = = 



B. Radiation-dominated era 

We now select Lagrangians that admit tracker solu- 
tions during radiation domination, w m = i. In order not 
to violate the nucleosynthesis bound, the energy density 
of fc-essence must be subdominant throughout the radia- 
tion era [3j, 



Ro > 0.99. 



(78) 



Admissible trackers may have a value Rq within the range 
0.99 < R < 1, or R = 1. A solution with < R < 1 is 
only possible with Lagrangians given by case 1, 

= i+y) lim ^ o( ^ = (79) 

We denote by v r the asymptotic value of v during the 
radiation era. Possible values of v r are determined from 
w,p(v r ) = and v r must satisfy 

C*(Vr) > |, £0(«r) ^ 0, U r j4 0. (80) 
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The corresponding value of Rq must respect the 
bound ([75]). 



Ro = l- 



9k 2 Q' : 



4e> 



> 0.99. 



(81) 



Solutions with Rq = 1 and possible in 

cases 3, 4, and 6. The hrst set of solutions is given by 

where 2 < a < ~. Admissible functions Ko((f>) decay or 
grow slower than any power of (f>, e.g. Ko(tf>) oc (ln<^)". 
Admissible values of v r are determined from the condi- 
tions 



2a 

■w<t>(v r ) = -y - 1, e^K) 7^ 0, u r ^ 0. (83) 
The second set of Lagrangians is 

Jf(« = ^, Um^)=0, lim^M = 0. 

(84) 

The possible values of u r are determined from w^ivr) = 
^, and the following conditions must be also satisfied, 

c 2 s (v r )>^, e$(v r )^0, v r ^0. (85) 

The third set of admissible Lagrangians is described by 
case 6 with n = 3, namely 

K U>) = ^5^. lim = 00, (86) 



09/4 ' 0— j-cxi 

Q(«) =Qiv 3 + o(v 3 ), Qi>0. 



(87) 



In this case, v r = 0. The solution of case 6 with n > 4 
cannot be used since the condition ([62]) cannot be satis- 
fied with w m — |. 



C. Dust-dominated era 

We now select the tracker solutions that exist for w m = 
0. In order to describe the late-time domination of k- 
essence, we must look for solutions with < —\ and 
i?o = 0. The possible trackers are cases 2, 8, and 10. 

In case 2, the Lagrangian must satisfy 



ff(*) = i±§M lim W )=0. 



We denote by Vd the asymptotic value of v during the 
dust era. The admissible values of Vd 7^ are determined 
from 



2w d 



(89) 



In addition, the following conditions must be satisfied: 
- 1 < w^(v d ) < 0, c 2 s {v d ) ± 0. (90) 
The second set of Lagrangians is for cases 8 and 10, 

K{4>) = ^M, lim KoM = oo. (91) 

This condition for K((j)) is satisfied, for example, by 
K((f) oc Q with a > —2. The value v d must be such 
that 

Q(v d ) < 0, Q> d ) = 0, (92) 
while we may have either vj, ^ or v d = 0. 

D. Viable scenarios 

Having listed all the Lagrangians that admit desired 
solutions in the radiation- and dust-dominated eras, it 
remains to determine the overlap between these classes 
of Lagrangians. By comparing the requirements on the 
functions K(<f>) and Q{v), we find only two possibilities 
for trackers in the radiation/dust era: case 1/case 2 and 
case 6 /case 8. 

The first set of Lagrangians (case 1/case 2) is 



K(<p) 



1 + K {(t>) 



lim K {(f>) 



0. 



(93) 



In the radiation era, the asymptotic value of v is given 
by v r 7^ such that 



w<f>(v r ) 



c 2 Jv r 



> 3, ^ 0, (94) 



and the dust attractor is given by Vd 7^ such that 
Eqs. (|89| - (|90|) hold. These Lagrangians describe the 
well-known scenario [2j where the fc-essence tracks ra- 
diation during the radiation era and eventually starts to 
dominate in the dust era. The function Q(v) must be 
chosen to satisfy the conditions of cases 1 and 2. Addi- 
tionally, one must exclude the possibility of a dust tracker 
(case 1, w rn = 0) by adjusting Q(v) such that the condi- 
tions of case 1 are not satisfied for w^(vq) = w m = @|. 

The second set of Lagrangians is described by case 
6/case 8. The function K ((f)) is of the form 



lim Ko((f>) = 00. 

<p — > 00 



The function Q(v) must be such that 

Q(v) = Q lV 3 + o(v 3 ), Qi>0. 



(95) 



Then the asymptotic values of v are v r = in the ra- 
diation era (where ~ i) and v d 7^ in the dust era 
(where w<p f» — 1). The value v d must be a root of Q'(v) 
such that 



Q(v d ) < 0, Q'{v d ) = 0. 



(96) 
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This scenario, however, has a fatal flaw. The attractor of 
case 8 requires that v > Vd at all times, while the attrac- 
tor of case 6 is realized at very small v « 0. Therefore, 
a transition from the first attractor to the second will 
necessarily involve values of v < Vd for which the theory 
is unstable since c 2 s {v) < 0. Hence, this scenario must be 
discarded. 

Thus we conclude that successful models of fc-essence 
are produced only by Lagrangians described by Eq. (|93)) 
under the conditions of case 1 and case 2. 



E. The existence of a superluminal epoch 

We have shown that the only viable fc-essence scenario 
is described by case 1 / case 2 of Sec. IV Al Now we demon- 
strate that in these scenarios c 2 (u*) > 1 for some value 
v * that is reached during the dust-dominated epoch. The 
argument is similar to that in Ref. [201 ] . 

Since Q{v r ) > and Q(vd) < 0, while Q(v ) is a mono- 
tonically growing function of v, we must have Vd < v r . 
In both scenarios of case 1 and case 2, the asymptotic 
fraction of the energy density Rq is equal to a certain 
function F of vq, 



R = F(v ) = 1 



9k 2 Q r 



4 vC 



(97) 



In case 2, F(vq) — due to Eq. (|4"5|) ; therefore, we may 
describe both cases 1 and 2 by a single function F(vq). 
We note that e^,(v) — vQ'(v) — Q(v) is a monotonically 
growing function of v because 



The equation of state parameter w^, (v) is a continuous 
function of v that satisfies 



= w^Vd) < 1 < Uty(«i). 



(102) 



Hence, there exists a value u* € [Kdj'fi] such that 
w^v*) > 1 and u^(v*) > 0. 

Finally, we show that c 2 (i;*) > 1 follows from the 
conditions u^i**) > 1 and w'^v*) > 0. According to 
Eq. d39j), we have 



Therefore 



(1 + W(jl ) (c 2 - W4) 



> 0. 



c 2 (v*) > W^,(v*) > 1. 



(103) 



(104) 



Since c 2 {v) is a continuous function, this demonstrates 
the existence of an interval of values of v within [w^iv] 
where c 2 {v) > 1. This superluminal epoch occurs during 
the dust-dominated era. 
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—sAv) = vC 
dv 



» > for v > 0. 



(98) 



Since e^,{v) > for every relevant value of v, it fol- 
lows that F(v) is a continuous function for these v. For 
a successful model of fc-essence, the radiation tracker 
must have F(v r ) > 0.99 and the dust tracker must have 
F(vd) = 0. During the evolution from the first tracker 
to the second, the value of v must traverse the interval 
[fd)W r ]. The condition F(vd) < F(v r ) implies (due to 
the continuity of F) that there exists a value v\ € [vd, v r ] 
such that F'(v±) is positive: 



F'(vi) = 



9k 2 Q'Q" 
2{vQ'-Qf 



> 0. 



(99) 



Since Q> > 0, Q" > 0, and e = vQ' - Q > for all 
v G [ud,fv], we can simplify this condition to 



vQ' 
2 



<0, 



or equivalently to 



> 1. 



(100) 



(101) 



Appendix A: ASYMPTOTICALLY STABLE 
SOLUTIONS 



The standard analysis of the /c-essence trackers 
(e.g. @) involves several simplifying but restrictive as- 
sumptions concerning the behavior of the solutions. A 
wider range of fc-essence models will be obtained if some 
of these assumptions are lifted. Let us therefore charac- 
terize the desired features of the cosmological evolution 
of fc-essence in a general manner. 

Scenarios of fc-essence are based on the assumption 
that the field <f> has an almost constant equation of state 
parameter (w<f,) during a cosmologically long epoch while 
another matter component dominates the energy density 
of the universe. Eventually, the fc-essence itself becomes 
dominant and plays the role of "dark energy," again with 
an approximately constant w$. It is important that the 
solution curves serve as attractors for all neighbor so- 
lutions. In that case, the value of at late times is 
essentially independent of the initial conditions. 

When the radiation-dominated epoch gives way to the 
epoch of dust domination, the behavior of fc-essence will 
change in a model-dependent way. However, it is techni- 
cally convenient to study the behavior of fc-essence under 
the assumption that the dominant matter component has 
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a fixed equation of state for all time. Then the existence 
of tracker solutions will be found by studying the asymp- 
totic behavior of the solutions at t — > oo (equivalently, at 
<t> — > oo). This is the approach taken in this paper. 

The evolution of fc-essence together with a single mat- 
ter component is described by the equations of motion 
(EOM) shown above as Eqs. in terms of the 

variables {v(<fr),e m {4>)}- We call a solution {«(</>), e m ((j))} 
asymptotically stable if Ws(4>) tends to a constant at 
(j) — ► oo and if all neighbor solutions (at least within 
a finite domain of attraction) also approach the same 
value of w,!,. In this section, we restrict our attention 
to asymptotically stable solutions with one matter com- 
ponent. Since reasonable values of are within the 
interval [—1,1], it is justifiable to ignore solutions where 
tends to infinity at late times. In principle, one could 
also have solutions where w c /,((j)) oscillates without reach- 
ing any limit as — > oo, but such solutions are of lit- 
tle physical interest since the value of w^, at the end of 
a given cosmological epoch will then be largely unpre- 
dictable. Since in the models under consideration is 
a function of v only, solutions v(cj>) that oscillate with- 
out reaching any limit are also excluded. Applicability 
of the effective field theory requires that the derivatives 
of 4> remain bounded; thus v = <f> cannot diverge to in- 
finity as <f> — > oo and must also tend to a constant value, 
v(4>) — > v < oo. 

There may also exist solutions with initially negligible 
but growing ratio e^/e m . Such solutions may have a 
stable behavior with an almost constant w$ for a finite 
(but very long) time, until the energy density of fc-essence 
starts to dominate. We do not consider such "transient 
attractors" in the present paper. 

Our main task is to deduce the possible fc-essence 
Lagrangians p(X, </>) that admit physically meaningful 
asymptotically stable solutions. We consider only La- 
grangians that have a factorized form (f!?0"]) . 3 We assume 
that the matter component has a constant equation of 
state parameter w m such that w m 7^—1. 

It will be convenient to use also the auxiliary variable 
R{4>) satisfying the EOM Since the values of R 

are limited to the interval [0, 1], any asymptotically sta- 
ble solution will necessarily approach a constant value, 
R(4>) — * Ro as (j) — > oo. The possible values of Rq and vo 
are yet to be determined; the cases when Rq or vq assume 
critical values (Rq = 0, Ro = 1, Vq = 0) will need to be 
treated separately. 

The general method of analysis is the following. We 
have a system of nonlinear EOM parameterized by a 
pair of functions K(<j>),Q(v); the general solution of the 
EOM is not available in closed form. Our purpose is 



3 Since the analysis uses only the properties of the Lagrangian 
in the asymptotic limit <j> — » oo, our results will apply to more 
general Lagrangians that have the form p = K(<f>)L(X) asymp- 
totically at large tj> and fixed X. 



to determine the functions K(<fi),Q(v) for which a so- 
lution of the EOM exists with the asymptotic stability 
property. We first assume the existence of an asymptot- 
ically stable solution {v(4>), e m (({>)} and derive the nec- 
essary conditions on the functions K(cj>) and Q(v) that 
admit such solutions (perhaps in more convenient vari- 
ables, such as {v ((/)), R(<fi)}). At this step, there will be 
many cases corresponding to different asymptotic behav- 
ior of v{4>) and R(4>). For instance, v(<j>) may tend either 
to a nonzero constant or to zero, etc. In each case, we 
then obtain the general solution of the EOM (with two 
integration constants) near the assumed stable solution 
(e.g. v{4>) = wo — A{<j)), with A{<j>) very small). At this 
point, it is possible to make simplifying assumptions be- 
cause we only consider the solutions in the asymptotic 
limit 4> —> oo and infinitesimally close to an assumed tra- 
jectory. We then investigate whether the general solution 
is attracted to the assumed stable solution. In this way, 
we either obtain sufficient conditions for the existence of 
a stable solution of an assumed type, or conclude that no 
stable solution exists in a given case. After enumerating 
all the cases, we will thus obtain necessary and sufficient 
conditions on K(<fi) and Q(v) for every possible type of 
stable tracking behavior. 

Let us begin by drawing some general consequences 
about the asymptotic behavior of stable solutions at <fi — > 
oo. Rewriting Eq. P^j) as 



de. 



-1/2 



((f)) 3n(l + w m ) 



d<p 2v(<j>)y/R{<t>Y 



(Al) 



and noting that the right-hand side of Eq. (|Aip is 
bounded away from zero, we conclude that £ m {4>) decays 
either as (f>~ 2 or faster at cf> — ► oo, depending on whether 
v(4>)y/R((f>) tends to zero at large <fi. In the following 
subsections, we consider all the possible cases. 

Based on the motivation for introducing fc-essence, we 
have assumed that w m ^ — 1. According to Eq. (jAljl . for 
w m < — 1 (phantom matter) the energy density e m will 
satisfy the differential inequality 



-1/2 



3/dl 



- w„ 



dt 



< -Ox, 



(A2) 



where C\ is a positive constant. Thus, e m (t) will reach in- 
finity in finite time regardless of the behavior of R(4>) and 
v(4>). However, this time can be quite long and the phan- 
tom behavior might be only a temporary phenomenon. 
Therefore, we will use the property w m 7^ —1 but avoid 
assuming that w m > — 1. 

In the analysis below, we will also use the following 
elementary facts: 

a) If a function F{<j>) monotonically goes to a constant 
at 4> — > 00, then F'(<f>) decays faster than </> -1 . This is 
easily established using the identity 



/>oo 

F(Q) - lim F{cf>) = F'{(j))d<j) 
4,^00 J 



< 00, 



(A3) 
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which means that F'((f>) is integrable at <fi — > oo. Hence, 
F'{<j>) decays faster than tfi -1 at <fi — * °°- 

b) If a function F(d>) is monotonic, then F'(d>) — > if 
and only if 



lim — -p- = 0. 

4>^rOO ffl 



(A4) 



This statement follows from the L'Hopital's rule in case 
F{<j)) — > oo, and is trivial in case i* 1 (0) has a finite limit 
at <p — ► oo. 



1. Energy density e m oc 2 and flo / 1, main case 

According to Eq. (|A1[) , the asymptotic behavior e m oc 
is possible only if v{<f))^ ' R(<f)) stays bounded away 
from zero as (f> — > oo, in other words if vq ^ and i?o ^ 0. 
We also assume Rq ^ 1, meaning that the energy density 
of fc-essence tracks the matter component; thus e^{<j)) oc 
d>~ 2 as well. It follows that H(<fi) oc y/e^,(<p) oc and 
then Eq. (JHI]) yields 



dR{4>) w Tt 
— r- — oc — 



w<t>{4>) 



(A5) 



Since i?(0) — > const, the derivative dR/d<p must decay 
faster than (j)^ 1 as d> — ► oo. Hence w^{4>) — ► w m as — > 
oo. This is the standard tracker behavior: the equation of 
state parameters of fc-essence and matter become almost 
equal at late times. 

By assumption, at large <f> the Lagrangian is factorized, 
p = K(<f))Q(v), and then we have 



w„ 



Q(vo) 



voQ'(vo) - Q(vo) 



(A6) 



This algebraic equation determines the possible values of 
vq for a given w m . (Tracker solutions of this type are 
impossible if this equation has no roots.) The property 

£<l>(<fi) oc <\T 2 becomes 



e (0) =K{ct>) (vQ'-Q) oc0- 2 . 
Generically one expects 

£<f>{vo) = vqQ'(vq) - Q(v ) ^ 0, 



(A7) 



(A8) 



and we temporarily make this additional assumption. 
Then we obtain 



K((f>) oc cf> 2 as (f> — > oo. 



(A9) 



This is somewhat more general than the function K(d>) = 
const • 4>~ 2 usually considered in fc-essence models. 

We may consider Lagrangians p — K(d>)Q(v) with the 
function K(d>) of the form 

= 1 + lim K (cf>) = 0. (A10) 



Let us now derive a sharp condition for the existence 
of an asymptotically stable solution {v((f>), R(<f>)} in this 
case. We use the ansatz 

v{<l>) = v -A{(t>), R{<f>) = R -B{4>), (All) 

where by assumption the unknown functions A(d>), B(<f>) 
tend to zero at d> — > oo. After deriving and solving the 
equations for A(d>) and B (</)), we will need to verify this 
assumption. 

Since the left-hand side of Eq. (|35[) is —A', it tends to 
zero faster than <f>~ . On the other hand, assuming that 
c 2 (vq) ^ 0, we find that the right-hand side of Eq. (f35|) 
contains leading terms of order (/> _1 , such as (In AT) ^ and 

\K. Hence, these terms must cancel, which entails 



3m 



£<t>( v o) 



2v Q 



1-Ro l + w<p(v ) 1 + 



2«o _ 264,^0) 
~ Q'(v ) ' 



(A12) 



Since vq is determined from Eq. (|A6|) . this condition fixes 
the value of Rq, 



Ro = l- 



9k, 2 Q' 2 (v q ) 
4 v Q'(v o )-Q(v o y 



(A13) 



The requirement that the values of i?o be between and 
1 further restricts the possible functions Q(v). Using 
Eq. (|A12[) , the condition R > can be expressed equiv- 
alently as 



< qTI^OT— 72- 

»« (l + w m ) 



(A14) 



No tracker solution is possible if this condition is violated. 

The equations for A(<fr) and B(<f) are now found by lin- 
earizing the equations (|35[) - (|36p . For brevity, we rewrite 
these equations as 



dv 

d\n4 
dR 
din d 



d\n < 



- A i( v )-7T-T ~ A ?(v)\ /— : 5^, (A15) 



1 - R 



-Aa(t>, R)y/l + K {<t>) (w m - w^v)) , (A16) 



where the auxiliary functions Ai, A2, A3 are defined by 

£</> 



A ( \ = c2 ^ V _ 

+ «Q"(«)' 

A 2 (v) = 3Kc 2 s (v)yJe^(v), 

3k 



(A17) 
(A18) 



A 3 (u, R) = —Ry/l-RJi^v) = ^ ^W' 



Note that Eq. (j A12|) is equivalent to 

A 2 M 



2Ai(«o) 



(A19) 



(A20) 
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Substituting the ansatz (|A11|) into Eqs. (|A15J> — (|A16[) . us- 
ing the identity (|A20[) , and keeping only the leading linear 
terms, we find 



dA 



(^+Jfo)A 1 (« )-a A- 
dine!. A 3i v o, Ro)w' 4> (v Q )A, 



dim 
dB 



AiHB 
1-iZo : 



where we defined the auxiliary constant ao by 



A 2 (u ) 



2Ai(v ) 



1 - Wn 



For the moment, we assume additionally that 



1 



i'o 



(A21) 
(A22) 

(A23) 



(A24) 

Differentiating Eq. (|A22[) with respect to ln0 and sub- 
stituting into Eq. (|A21|) , we find a closed second-order 
equation for A(<p), 



d 2 B 



dB 

aa-p. — T + PqB = 7o 



dK Q 
din (A 



d(ln0) z din 
where the constant coefficients Poilo are defined by 
7o = Ai(v )A 3 (v , Ro)w' <j> (v ) 



, (A25) 



= 2- 



1 + 



-w m R (1 - 5 ) 



A) 



7o n c s {v )-w m 2 
1 fT = — TT W m^0- 

1 - -Kn 1 + Wm 



(A26) 
(A27) 



The general solution of Eq. (|A25|) is the sum of an 
inhomogeneous solution and the general solution of the 
homogeneous equation. Homogeneous solutions are sta- 
ble if both roots Ai 2 of the characteristic equation 



A + a Q X + f3 = 
have negative real parts, 

Re(Ai)<0, Re(A 2 )<0. 
This will be the case if 

a > 0, A) > 0, 
which is equivalent to the conditions 

\w m \ < 1, C 2 s (v a ) > W m . 



(A28) 
(A29) 
(A30) 
(A31) 



An inhomogeneous solution of Eq. (|A25[> can be ex- 
pressed as 



B{cj>) = Bi 
5 2 (0) 



B- 



(A32) 



To 



Ai — A2 

7o 
A? — Ai 



0- Al - 1 (0i^ + if o )d0, (A33) 
oT^" 1 (cj)K' Q + K ) d0. (A34) 



Since the function -Ko(0) tends to zero at — > 00 by 
assumption, the inhomogeneous solution also tends to 
zero at — > 00 as long as the condition (|A29| holds. 
This is straightforward to show by assuming an upper 
bound 



\<f>K' Q + K Q \<M for all > M , 



(A35) 



where m can be chosen for any M > 0. Then the inho- 
mogeneous solution 5(0) is bounded for > 0m by 



1 5(0) I < const • M + const ■ 



const ■ 



(A36) 



which means that 5(0) — > at — > 00. 

Under the same assumptions, the function A(<f>) will 
have the same behavior at — ► 00. We conclude that 
asymptotically stable solutions {w(0), 5(0)} approaching 
{vo,Rq} exist under the assumption c^(vq) 7^ and the 
further conditions flXB]), CM)), (|AT2]l . JX24}, and (fX3T|) . 4 
These conditions are similar to those derived in Ref. Q 
under a more restrictive assumption K (0) = const • 0~ 2 . 
Let us now investigate whether these assumptions can be 
relaxed further. 



2. Energy density e m oc and i?o 7^ 1, marginal 
cases 

The last assumption used in the derivation of the sta- 
bility condition (|A31|) was Eq. (|A24j) . If c 2 (v ) = w m 
while all the other assumptions hold, we have w'Avq) = 
and the equation (|A22[) for 5(0) is modified. We may 
then rewrite Eqs. (|X2T]) - ((A22|) as 



dA 
dine' 

dB 
d\nd 



(0X o + X o )A 1 (w o ) -a A 
0{A 2 ). 



AiM5 
1 -Ro ' 



(A37) 
(A38) 



Differentiating the first equation with respect to In 0, we 
obtain 



d 2 A 



d(ln0) 



dA 

, - 0(0if o )" Ai(w ) - a a — — 



dln< 



(9(A 2 ). (A39) 



The second-order terms 0(A 2 ) can be disregarded for the 
stability analysis. Since the characteristic equation 



A + a X = 



(A40) 



has a zero root, the general solution {A(<fi), 5(0)} will not 
tend to zero at — > 00. Hence, no asymptotically stable 
solutions exist when the condition (|A3ip is violated. 

Another assumption, c 2 (vq) 7^ 0, was used to derive 
Eq. (|A12p that determines the allowed value of 5o- Let 



4 This is case 1 in Sec. |VX1 
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us briefly consider the possibility c 2 (v ) = 0. (We note 
that v 7^ Vq on actual trajectories, so stability will hold 
as long as the trajectories v(4>) do not reach the regime 
e» < 0.) If 



Q'jvo) 
v Q"(vo) 



0. 



then Q'(vq) = and the asymptotic equation of state is 



«ty(«o) 



Q(v ) 



voQ'(v ) - Q(vo) 



(A42) 



as long as Q(vq) ^ 0. However, we assumed a matter 
component with w m =/= — 1, and so we discard the possi- 
bility that Q(v ) ^ 0. If, on the other hand, Q(vq) = 0, 
then we must also have e<p(v Q ) — 0. Thus c 2 (v Q ) ^ is 
justified given that e^vq) ^ 0. 

Relaxing the assumption e^{vq) ^ requires some 
more work. If es{vq) — 0, then we cannot conclude 
that K(<p) oc <\T 2 at <f> — ► 00; the function K{<f>) re- 
mains undetermined even though we know that s^(<p) = 
K((j))e < f,(v) oc </>~ 2 . The analysis after Eq. (|A6[) needs to 
be modified as follows. The finiteness of w^,, 



w<b\vo) = hm _ < 00, 



(A43) 



requires that Q(fo) = and thus (since vq ^ 0) also 
Q'(vo) = 0. In general, we may suppose that Q(v) has 
an expansion 

Q(v) = ^(v- v Q ) n [1 + 0(v - v )\ , (A44) 

where Qo is a nonzero constant and n > 2. In this case 
we have the expansions 

= Qo (v - vof- 1 [1 + 0(v - «o)] , (A45) 
w - w 



nv 
« - wq 
(n - 1) uo 



[1 + 0(v - v )\ , 



[1 + 0(v - «q)] 



(A46) 
(A47) 



It follows that W0(«o) = 0, so the only possibility for 
tracking is w m = 0. Also, the only admissible solutions 
are those with v(cf>) > Vq, meaning that A(<f>) < and 
Qo > 0. Let us now perform a stability analysis of these 
solutions. Substituting the ansatz (|A11|) into Eqs. (|3"d])- 
(|36p and keeping only the leading terms in the perturba- 
tion variables A{<f>) and B ((/)), we obtain 



dA {-A) K' 



3 k 



n-XK (n-l)uoV 1 ~ R ° 



(n+l)/2 



(A48) 

^ = -^i?ov / T^v^(^(-A) ( " +1)/2 . (A49) 
Wy) ^^0 

For the purposes of a stability analysis, it is sufficient to 
note that Eq. (|A48[) does not involve B((f>). One can solve 



Eq. (IA48P explicitly for A(<p) and find such K(4>) that 
the general solution for A(<fi) tends to zero at <f> — > 00; for 
instance, K(<p) oc <jf with r > —2. However, the general 
solution for B((f>) is 



(A41) B(<i>) = B - const • / y/K($) (_A) (n+1)/2 d(f>, (A50) 



where Bo is an arbitrary integration constant. It fol- 
lows that B(<fi) will either diverge or tend to an arbitrary 
constant of integration at <p — > 00. Hence, the general 
perturbation will not tend to zero at large <f>. We con- 
clude that no asymptotically stable solutions exist when 
£<t>{vv) = 0. 



3. Energy density e m oc <j> 2 and Ro = 1, main case 

We use the ansatz R{4>) = 1 — B(<f>), where the func- 
tion B(4>) is positive and tends to zero monotonically as 
(j) — > 00. Since dR/d(p > 0, it follows from Eq. (|3"T|) that 
w m < w,p(v((j))) for all sufficiently large 0. Thus, any 
asymptotically stable solutions will necessarily satisfy the 
condition 



W m < Wj,(v ). 



(A51) 



Since R — ► 1 as <p — > 00, we have £tot(<^) oc £ m {4>) oc 

(,-2 



', so we may write 

£tot(0) ~ So^ _2 , 



00, 



(A52) 



where i?o is a nonzero constant. The value of Eq can be 
related to other parameters by using Eq. (|25p . rewritten 
as 

rf In e m 



3 k 



dlncf) v 
which yields, in the limit 

2 = - V 

Expressing £ to t through £0, we have 

£0 _ e<p(v)K(<p) _ 



— y/(/) 2 e m (l + 
t> — > 00, 
S (1 + w m ) . 



So 1 ^ 2 ~ £tot 



I - i? 



hence 



B(0 



£ (t;)0 2 g(0) 
So 



(A53) 



(A54) 



(A55) 



(A56) 



We now assume that i^(vo) ^ 0; the case £0(^0) = 
will be considered later. If £0(^0) 7^ 0, it follows that 



Rewriting Eq. ([31]) as 

dln(l -R) _ 3kR 

dhxd) v 



(A57) 



V 4> 2 £tot {w m - w<j,(v)) (A58) 
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and substituting Eqs. (|A52p and (|A57[) . we find for large 
d\u{4> 2 K((j))) 3KA/Eb w m -w$(v) 



(w m - w<p(v)) = 2- 



d\n(f> vo 1 

(A59) 

It is now clear that the possible asymptotic behavior of 
K((f>) at <f) — ► oo depends on whether W^,[v) tends to w m 
at large 4>, i.e. on whether or not w^fao) = w m . (We 
note that the value of vq is yet to be determined by the 
analysis that follows.) 

Considering the interesting case w^(vq) ^ w m , we find 
that the right-hand side of Eq. (|A59[) tends to a negative 
constant as <j) — > oo. Denoting that constant by —fj,, 
where 



>0, 



(A60) 



and integrating Eq. (|A59p , we infer the following asymp- 
totic behavior of K(<p), 

K(<P) cx (p-z-^Ko^), ^ oo, (A61) 

where -Ko(0) is an auxiliary function that satisfies 

dhxK (<p) 



lim 

d> — >oo 



dhi ( 



This condition is equivalent to 



lim 



ln< 



0. 



(A62) 



(A63) 



Thus, the function Ko((f>) may go to a constant at large 
4>, or may grow or decay slower than any power of 0; 
examples of admissible functions Ko((j)) are 

K (4>) = (ln<f>) p ; Xo^HexptdQn^), |*| < 1. 

(A64) 

With any such Ko((fi), solutions of the currently con- 
sidered type are possible only for Lagrangians p = 
K((j>)Q(v) with 



where 



K{4>) = r 2a Kv(4>), 

2 + /it _ 1 + w^{vq) 
2 1 + w m 



> 1. 



(A65) 



(A66) 



For a given Lagrangian of this type, the possible values 
of v are fixed by Eq. (|A66|) . If Eq. (|A66[) is not satisfied 
for any such vq, solutions of this type do not exist. The 
value Wj,(va) is determined by Eq. (|A66|) as 



w<t>(v ) = (1 + w m ) a — I. 



(A67) 



Since w m ^ —1, we must have w$(vo) —1 also. 

It remains to investigate the asymptotic stability of the 
general solution. Since #(</>) must satisfy Eq. (|A56|) . we 
may write an ansatz 



B{4>) 



E 



(A68) 



where C((j>) is a new perturbation variable. Hence, we 
substitute Eq. (|A65[) together with the ansatz 



v(<fi)=v -A(4>), 
R{<t>) = 1 - 

£tot(«, 4>) = E (, 



Eq 

_ 2 1 



E^vq) 1 + C{4>) ' 



(A69) 
(A70) 

(A71) 



i nto E qs. (|35| and JMSI- Using Eqs. (TA"54"|) , !|A"60|) , and 
(|A65|) . we obtain at an intermediate step the equations 



dA 



2a , . / 

-— + (bxK y 

9 

(j)- 1 2v 



Ai(») 
A 2 («) 



1 dC 
1 + C ~d4> 



VTTci + 

^f 1 -(hxK )' 
H^- 1 A 4 (u) 



(A72) 



v / r+CA 4 (t> ) 



[1 + 0(0-")], (A73) 



where the functions Ai(u) and h.2{v) were defined by 
Eqs. (|A17[) - (IA18[) , while the new auxiliary function h.±(v) 
is defined by 



A 4 (w) = Sn^Je^v) 

In the present case, the identity 

2v 



~ w 4>( v ) 



2aA 1 (v ) 



A 2 (v ) 



1 



w m Znyje^va) 



(A74) 



(A75) 



holds due to Eq. (|A"66]l . 

We now linearize Eqs. (|A72|) ~ (|A73|) with respect to the 
perturbation variables A and C . To simplify the lin- 
earized equations, we use Eqs. ([3"9"| . (|A54|) . and the defi- 
nition (|A66[) of a. (We note that Ai(i>o) ^ 0; otherwise, 
we would have c 2 s {vq) = 0, which contradicts the earlier 
assumptions e^,(uo) ^ and w<f,(v ) ^ — 1.) After some 
algebra, we find (to the leading order) 



dA 
dlnd 



2aA 2 (v ) 



AiM 
A 2 (v) 



A 





+ Al(l! ) 


~d\nK Q 
dlncp 


-aC 




dC 
dhi(j) 


d\nK 
dhicp 


+ \*c- 


KM 
hM A 4 M 



(A76) 



A 



s<i>(vo) 



(A77) 



This is an inhomogeneous linear system for A(4>) and 
C(<f). The analysis of the asymptotic stability is similar 
to that after Eq. (|A29[) . Since all the inhomogeneous 
terms are decaying at <j> — > oo, it suffices to require that 
both the eigenvalues of the homogeneous system have 
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negative real parts. For a homogeneous system of the 
form 



dA 

d\n4 
dC 
din c/ 



: Pi A + foC, 

7iA + 7 2 C, 



the characteristic equation is 

A 2 - (fa + 72 ) A + (/3 X 72 - /3 2 7i) = 0, 
and the stability conditions are 

Pi + 72 < 0, ^172 - P2J1 > 0. 



(A78) 
(A79) 

(A80) 
(A81) 



Presently, the constants pi , P2, 71 , 72 can be read off from 
Eqs. (|A76jl - ljA77|) ; simplifying, we obtain 



Pi 



7i 



1 -w^jvo) 
-a- ; — r . Pi 



a c 2 {v )vo 
](v ) - W^Voj) (1 + w m ) 



(A82) 



v oC 2 s {v ) 



1 - w,p(v ) 
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(A83) 



The stability conditions (|A81|) can be simplified to 

c 2 s (v ) - W^Vq 



1 - W^Vp) jj_ 
1 + w^(vq) 2a 



Wc/>(vo) - W n 



> 0. 



(A84) 



Since w,p(vo) > w rn for solutions of the present type, 
while 2a = /1 + 2, the stability conditions (together with 
the condition u^(vo) > w m ) are 

w m < w^vq) < ~r—~ ' c l ( v o) > w^vq). (A85) 

_L ~\~ fJj 



Using Eq. (|A67[) . we can transform the first of these con- 
ditions into a condition for a: 



1 < a < 



1 



1 



1 + Wr, 



c 2 (v ) > W^Vq). 



(A86) 



The first inequality above will define a nonempty interval 
of a only if \w m \ < 1. These are the final conditions for 
the asymptotic stability of the solutions obtained under 
the assumptions s,f,(vo) 7^ 0, w^,(v a ) 7^ w m , and (|A66|) . 5 



4. Energy density e m cx </> 2 and Rn = 1, marginal 



The analysis in the previous section used the assump- 
tions £^(^0) 7^ and w^Vq) 7^ w m . In this section we 
lift these assumption, in the reverse order used. 



This is case 3 in Sec. IV Al 



If w$(vo) = w m while £^(^0) 7^ 0, then we may con- 
tinue the arguments starting with Eq. (|A59[) . Note that 
Eqs. (|A54jl and (|A57|I still hold. Since the right-hand 
side of Eq. (|A59|) tends to zero at </> — > 00, it follows that 



hm — — = 0. 

0-^00 dlnffl 



This condition is equivalent to 



hm = _ 2 . 

rA— >oo In q> 



(A87) 



(A88) 



Also, according to Eq. (|A57[) we can have B(<fr) — > only 
if 



lim 4> 2 K(4>) = 0. 



(A89) 



So the function if (</>) cannot have a power-law asymp- 
totic other than </>~ 2 ; more precisely, for any e > and 
for large enough </> we must have 



K(<f>) < <fr 



K(4>) > 4> 



-2-e 



(A90) 



However, a non-power law asymptotic behavior at <p — > 
00 is still admissible, for instance K{<j>) oc </> -2 (ln</>) s , 
where s > to allow B(4>) — ► according to Eq. (|A57|) . 
Rather than assume a particular form of K{4>), we 
will perform the analysis for arbitrary K(cj>) satisfying 
Eq. Mil . 

We again use the ansatz (| A68I) to linearize Eqs. ([55]) 
and (|A58|) . After some algebra, we find (to the leading 
order) 



dA 
dlnd 



dC 
(fine/ 



din 1 



din ( 



Ai(v) + 



A 2 (u) 



E 



VTTcy e^(«o)' 
(A91) 



dln^K) 
d\nd> 



A 4 (u) 



(A92) 



where the auxiliary functions Ai(v), A 2 (w), and A^v) 
were defined above by Eqs. (fATTf . (fAT8)) . and [fATIj) . 
Since ii>0(vo) = w m and c 2 (vq) 7^ 0, we have the relation- 
ship, 



2Ai(v ) = A 2 (u )i 



En 



£4>( v o) 



7^0, 



(A93) 



and then, assuming for the moment that A'^vq) 7^ 0, we 
can linearize Eqs. (|A91j) - (|A92j) as 



dA -2A 2 M^ 1(W) 



din 



dC 
dine/ 



A a («) 



din (e/. 2 ^) 

din 
2u A^(w ) 



Ai(wo), 



3k v /£0(w O ) 



dlnfc/) 2 ^ 



din < 



(A94) 
(A95) 
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The stability analysis proceeds as before, since all the 
inhomogeneous terms are decaying at </> — > oo. The re- 
sulting conditions are simplified to 



> 0, 



>0, 



■ w ri 



and further to 



< 1, cf > w ri 



(A96) 



(A97) 



These conditions are the same as the standard stability 
conditions for a tracker solution. Under these conditions, 
a tracker solution with w^,(vq) = w m exists as long as 
K{4>) satisfies Eqs. (jM8l) - (jM9]) . 6 

Finally, we analyze the case £<f,(vo) = 0. Since 
Eqs. (|A52j) . (|A54|) . and (|A56|) still hold for an asymptot- 
ically stable solution, we are motivated to use the ansatz 



R(4>) 



va-A((f>), 
1 - B(v,</> 



B(v, 



_ i<t>(v) ,, 



E 



^K(<P) (1 + C{<f>)) 



(A98) 
(A99) 

(A100) 



We first derive the exact equations of motion for the vari- 
ables A(<f>), C(4>) from Eqs. §5$ and $MB'- 



dA 
dlncp 

riln(l + C) 
din 6 



vc 2 s (v) dlnK 2vq cl(v) 



w^v) dlncj) 
2«o 



1 + w m Vl + C 

(A101) 



-A , . _ 

v v \ vl + C 

2v B(v,4>) w r , 



1 



w<t,(v) 



(A102) 



Then the stability analysis consists of checking that the 
general solution involves functions A(cj>), B(v, cf>), C(cj>) 
that decay as cj> — > oo. Since vo ^ 0, the expan- 
sions (|A45I) - (|A47|> hold with n > 2; we note that A < 
to guarantee c 2 > 0, and that Wj,(vo) = 0. The leading- 
order terms in Eq. (|A101|) are 



d\A\ 



\A\ fdlnK(cj>) 



dlncp n — 1 \ dlncj) 
and the general solution is 

\A\ = A [cj)^K(q)) 
Since n > 2, solutions A(cj>) decay at 
K{4>)cj) T +^ — > oo, <f> — ► o 



(A103) 



(A104) 
oo as long as 
(A105) 



The function B(v,cj>) is then expressed as 

B{v,4>) = ^—if>^(l + C{il>)), (A106) 

and its decay at cf> — > oo requires that — 1 < w m < 0. 
The leading terms of Eq. (|A102|) are 



dC _ 2_ 
dlncj) vq 



A-C 



l-\c] B(v,<fi)-^- 
I 1 + w„ 



Since the homogeneous solution C(cj>) oc cj)^ 1 decays as 
a power of cj>, while the inhomogeneous terms all decay 
at cj> — > oo, the general solution C(cj>) will also decay at 
cj> — ► oo. Thus, we find a family of asymptotically stable 
solutions corresponding to a value v such that Eq. (|A44[) 
holds, in case w m < and for Lagrangians with K(cf>) 
that either does not decay at large cj>, or decays slower 
than<r 2/(1+,Um) . 7 



5. Domination by fc-essence, vq ^ 0, main case 

We now consider the case Ro = 0. In this case, the 
matter component becomes subdominant at late times, 
so £tot ~ at cp — > oo. According to Eq. (|2"6"|) . we have 
at late times 

4r £ <p( < i ) ) = -— v^tot (i + «v)£0 



3k 3/2 M . , » 



(A107) 



thus the asymptotic behavior of £0(0) depends on 
whether or not w^vq) = — 1. With «o 7^ 0, one can 
have w,p(vq) = —1 only if Q'(vq) = 0, which entails 



lim 



lim 



0. 



vo v ^ v o Q"{v) vo v^v (lnQ'(u))' 

(A108) 

Let us postpone the consideration of the case c s (vq) = 
0; thus, presently we have w^vq) 7^ — 1. In that case, 
the asymptotic behavior of E<p(cj>) and etot(^) can be ex- 
pressed as 



etotO) ~ £</>(» ~ E cj) 2 , 
where the constant Eq is given by 

2v 







1 + 'w ( j > (vo) , 



(A109) 



(A110) 



due to Eq. (|A1Q7|> . We use Eqs. (|35] ) -([36 ]l to describe 
asymptotically stable solutions. Since on such solutions 
R{cj>) approaches zero while remaining positive, we must 



This is case 4 in Sec. IV Al 



This is case 5 in Sec. IV Al 
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have w<j,(v) < w m at late times. Computing the limit of 
Eq. ([55]) as (j> — > oo and using Eq. (|A1 10|) . we find 



= lim 



dv 



fi—>oa din ( 



= lim 4*c s (v) 

(j) — >oo 



1 + W^(v ) 0->oo 



lim Cg(v) 



1 + Wt/,(v) 

d\nK{4>) | 2 
din 6 



(Alll) 



The right-hand side of Eq. (|A111[) can vanish at (f> — > oo 
if, for instance, c 2 s (vo) = 0. We postpone the consider- 
ation of the case Cg(vo) = and presently assume that 
c s( u o) 7^ 0, which (together with w ^ 0) also implies 
£<j>{vo) 7^ 0. Then e^{<f>) oc (f>~ 2 entails K{<j>) oc (f>~ 2 at 



</> — » oo; accordingly, the right-hand side of Eq. (|A111|) 
vanishes at (/> — ► oo due to 



lim — = -2. 

<4— >oo dmd> 



(A112) 



By absorbing a constant into Q(v) if necessary, we may 
express K(<f>) as 



K{4>) = 



lim X- o (0) 

0^-00 



0. 



(A113) 



This is the familiar form of the function K(cj>), shown by 
Eq. (fATOl) in Sec. [ATI 

For these K(<p), the condition (|A110[) becomes 



3kv £Vi 



= ZK^Ji^vo) = — 



2wo 



^(/Uo) ' 



(A114) 



which is an equation for determining the admissible val- 
ues of vq. For these vo, we linearize Eqs. (f3"5| ~ (f3l)|) using 
the ansatz 



v a -A{cj>), R = B{cj>), 



(A115) 



where A(<p), B(<f>) tend to zero as <f> — > oo. The manipula- 
tions with Eq. ([33)) are the same as those in Sec. IA 1[ the 
result of the linearization is quite similar to Eq. (|A21[) 
with i?o = and without the relationship w^(vq) = w m , 



dA 



(cj>K +K )A 1 (v ) 



din 

The linearized form of Eq 

dB _ 3k^/Eo~ 
d In <fi vq 



1 - W^Vp) 

1 + w<p(v )' 



A 



A 1 (v )B. 
(A116) 



is 



{w m - W<fi(Vo)) B 



_ 2 Wm ~ W<t,{vo) 
1 + W^Vq) 



(A117) 



Since the equation for B(<fi) does not involve A(<f)), and 
since w m > w<f,(va), all solutions B((f>) decay, and thus 
all solutions A{<j>) also decay as long as 



1 - w<f,(vo) 
1 + w,f,(v ) 



>0, 



(A118) 



which is equivalent to |u>0(vo)| < 1- Therefore, solutions 
are asymptotically stable under the conditions (|A1 14|) . 
|t«0(t>o)| <l,w m > w^vo), and c s (v ) ^ 0. 8 



6. Domination by fc-essence, vq ^ 0, marginal cases 

In this section we continue considering the case Rq = 0, 
Vq 7^ 0, and examine the possibility that c 2 (vo) = 0. In 
that case, we have Q'(vq) — as well, which fixes ad- 
missible values of Vq. There are two further possibilities: 
either Q(vo) ^ or Q(v ) = 0. 

If Q(vq) = Qo 0, then Q(v) can be expanded about 
V — vo as 

Q(v) =Qo + Qi{v- voT [1 + 0(v - v Q )} , (A119) 
where n > 2. One readily obtains the expansions 

= -Qo + nvoQi (v - wo)" _1 [1 + 0(v - v )] , 

(A120) 

/ \ 1 , nvoQi , ,n-l h v, 

w^(v) = -H — — (w-wo) [1 + 0(u - uojj , 



-Qo 

1 V — Vq 



DO ""I 



[1 + 0(v - v )} . 



(A121) 
(A122) 



It is clear that one must have K(4>)Q < due to the 
positivity of the energy density. For convenience, let us 
assume that K{<f>) > and Qo < 0. Substituting the 
expansions above into Eqs. (|33|) - ([5rj|) together with the 
ansatz (| Al 15|) and neglecting the subleading terms, we 
obtain 



dA 



dB 



1 



-A 



Vq 1% — 1 



— — — — T + 3«VIQo I , 
nQi {-A) 

(A123) 



3k 



— = B^\Qo\K(4>) (w m + 1) . 

a<p vo 



(A124) 



Since these equations are uncoupled in the leading or- 
der, the stability analysis is performed for each equation 
separately. Integrating Eq. (|A124[) , we find the general 
solution 



B(<j>) = exp 



Co 



vo 



{Wm + 1) / yjK{4>)dcp 



(A125) 

where Co is an integration constant. The general solution 
B(4>) will tend to zero if and only if J J K{(j>)d(j> diverges 
as 4> — > 00 and w m > — 1. Let us temporarily denote 



x(<P) 



X — > oo as <fi — > oo. (A126) 



8 This is case 2 in Sec. |VX1 



18 



Then we rewrite the first equation as 



dx 



\Qo\ K>(4>) 



nv Qi K 3 / 2 

(A127) 

(Note that we must have A < on solutions, due to 
the requirement of positivity of c 2 .) The general solution 
A(x) can now be written explicitly, but it suffices to ob- 
serve that A(x) will approach zero as x ~~ * 00 if an d only 
if 



lim 



-2 lim 



-IT 



-1/2 



0. 



This condition is equivalent to 



lim <j>y/K{<t>) = OO. 

0— ^OG 



(A128) 



(A129) 



Note that the condition (j A126|) follows from that of 
Eq. (|A129p . To verify this more formally, consider a func- 
tion K{4>) such that j°° y/W^)d(t> < 00. Then K 1 ' 2 ^) 
necessarily decays faster than </> _1 at <p — > 00, and so 
K ~ x l 2 grows faster than (f> at <fi — > 00. Such if(0) cannot 
satisfy Eq. (j A129|) . Therefore it is sufficient to impose 
only the condition (|A129j) . This condition is satisfied, 
for instance, by functions K{(f) oc <j) s with s > —2. Thus, 
we conclude that the solution with Rq = is asymptoti- 
cally stable under the condition (|A129p and assumptions 
Q(vo) ^ 0, Q'(v ) = 0. 9 

It remains to consider the case R = 0, Q(v ) = 
Q'{ v o) = 0. In that case, similarly to that discussed in 
Sec. lA2[ we may use the expansions (|A44I) - (|A47|) . It fol- 
lows that ^(^o) = 0. With the ansatz v((f>) — v — A(cj)), 
we find that A{<p) < on physically reasonable solutions. 
Then the leading terms of Eq. (|35p are 



d(-A) (-A)K' 3ny/Qo 



n — 1 K (n — 1) «o 



VW)(-A) 



(n+l)/2 



(A130) 



Since this equation is independent of B, it suffices to en- 
sure that A{(f) — ► as (j) — > 00 and subsequently consider 
the general solution for R(4>). The general solution for 
A((j>) can be easily found by rewriting Eq. (|A130|) as 



(-A) 



d 

d~4> 

We find 

f4)(»-l)/2 



-(n-l)/2 K -l/2 



2v 



3k\/Qo 
2v 



3k\/Qo 0-00 ^/K{4>) ' 



(A131) 



(A132) 



where (f>Q is a constant of integration. It follows that 
A{(j)) — > as 4> —* 00 if K{4>) is suc h that c/PK^) — > 00. 
Under this assumption, we find that 



£0 = K^e^v) oc 4> 2 , — > 00, 



(A133) 



as it should according to Eq. (|A109|) . Now we analyze 
the general solution for R((f>) Then the leading terms of 
Eq. ([36]) are 



dR 
dhid 



= -R- 



A \ ( A 
m + = -2R [w m + 

(A134) 

where we used Eq. (IA1 10|) . If w m = 0, the right-hand side 
above is always positive and (since R is always positive) 
the general solution for R{<fi) cannot approach zero. If 
w m 0, the general solution for R(4>) is 

R(4>) oc cf>~ 2Wm as (j) — ► 00. (A135) 

It follows that the general solution R((f>) — ► at <j) — * 00 
as long as w m > 0. We conclude that an asymptotically 
stable solution exists in case Q(vq) — Q'(vq) — if w m > 
and 4> 2 K(<j)) — > 00 as (f> — > 00. The admissible functions 
K((f) are, for instance, K((j>) oc 4> s with s > — 2 



10 



7. Slow motion (vq — 0), main case (Q(0) 7^ 0) 

Previously we have been assuming that vq 7^ 0. Now 
we turn to the case vq = 0, which means that the velocity 
(f> = v(4>) of the field </> approaches zero, albeit sufficiently 
slowly so that 4> still reaches arbitrarily large values at 
late times. We will now obtain the conditions for the 
existence of asymptotically stable solutions with v(<p) — > 
at (j) — > 00. 

The finiteness of the speed of sound at v — > 0, 



lim c\ 

D->0 * 



(v) = lim 



Q'(y) 

vQ"(v) 



< 00, 



(A136) 



requires that Q'(0) = 0. Since the important quantity 
£<t>(v) = vQ 1 — Q approaches — Q(0) at late times, it is 
useful to distinguish two possibilities: Q(0) ^ and (less 
generically) Q(0) = 0. In this section we consider the 
generic case, Q(0) = —Qo ^ 0. Positivity of the energy 
density requires that K(4>)Qq > 0, and we will choose 
K{4>) > and Qo > 0. 

Under these assumptions, we may expand the function 
Q(v) near v = as 



Q(v) = -Qo + Qiv" [1 + 0{v)} , 



(A137) 



where n > 2 is the lowest order of the nonvanishing 
derivative of Q(v) at v = 0, and Q\ > because Q(v) is 
a convex and monotonically growing function of v. Other 
relevant quantities are then expanded as 

e*(v) = Qo + (n - 1) Qiv n [1 + 0(v)\ , (A138) 
nQi 



w4v) = -l + -^v n [l + 0(v)}, 
wo 



C 2 s(v) = 



n - 1 



[1 + 0(v)} 



(A139) 
(A140) 



This is case 8 in Sec. IV Al 



This is case 9 in Sec. IV Al 
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It follows that the only possible equation of state is 
«u(0) = — 1, indicating a possible de Sitter tracker so- 
lution. 

The equations of motion (|3"5)) - (j3l)|) become (neglecting 
terms of order v) 



dv 
d0 

dR 



1 



Qo K' 



nQiV 7 ' 



K 



3k 



K(4>)Q 
1-R 



3k 



Ry/T^Ry/ K{<j>)Q {w n 



1 



(A141) 
(A142) 



The first step is to investigate the possibility that R{(j>) — > 
1 at large <f> (we will hnd that this possibility cannot be 
realized). We note that for w m > —1, the right-hand 
side of Eq. (|A142[) always remains negative. Thus, for 
w m > —1 the general solution R(<f>) cannot tend to 1 at 
4> — ► oo, regardless of the behavior of K(<j>) and v(<p). In 
case w m < — 1, we need to do more work to establish that 
there are no asymptotically stable solutions with Rq = 1. 

Substituting the ansatz R(4>) = 1 — B(<f>) into 
Eq. (|A142[) and assuming that B — > 0, we obtain (omit- 
ting terms of order v and B) 



dy/B 



3 k 



= - — yjK(<t>)Q \l + W m \. 
Iv 



Changing the variable from (f> to \ defined by 



x(4>) 



we find 



dy/B 
dx 



(A143) 



(A144) 



(A145) 



There are now two possibilities: either the integral in 
Eq. (|A144[) diverges at 4> — > oo, or it converges. Accord- 
ingly, either x ~ > 00 or X ~ v Xo < oo at <f> — ► oo. In case 
X — > oo at <p — ► oo, we would have 



d^fB 

x^oo dx 



lim 



0. 



Since the right-hand side in Eq. (|A145|I tends to infinity 
at — > oo, the case x ~ * 00 is impossible. Thus, the 
integral in Eq. (IA144[) must converge at — > oo. It 
follows that K(<fi) — > faster than (j)~ 2 at — > oo, and 
then we may express K((f>) through an auxiliary function 
K (4>) as 



K{4>) = r 2 K {ct>), 



lim K Q {4>) = 0. 

(/> — too 



(A147) 



Further, we rewrite Eq. (j A141|) as 



dv 
din c 



Qo 



nQiv r 



2- 



dlnATo 
dhicb 



-3k 



W)Qo 



(A148) 



By construction, 



lim 12-^° 

<#>-+oo V am 



> 2 



(A149) 



(the limit might even be positive infinite if Kq tends to 
zero sufficiently quickly). Hence, under the assumptions 
v((f>) — > and B(cj)) — > we must have 



Qo 



lim 

4>^oo nQiv n 



2 - 



dlnK 
din (A 



(A150) 



It then follows by taking the limit <f> — > oo of Eq. (|A148[) 
that the two terms in the brackets must cancel while both 
approach infinity. Therefore, at large <j) we must have the 
approximate relationship 



K<t>) 



i 



y/B{fj SnnQi ^K Q {<j>) \ dln^ 

(A151) 

Due to Eq. (|A149[) . the auxiliary function M(cf>) defined 
by Eq. (|A151|) has the properties 



lim M^y/KoJf) > 



lim M{4>) = +00. 

0— >00 



(A152) 

(The first limit may be positive infinite.) Using the func- 
tion M(<f>), we may rewrite Eq. (|A148|) as 



dv _ 3Ky/ K Q ((j))Qo 



d\ii( 



1 



M 



(A153) 



Expressing \fB through v using Eq. (|A151j) and substi- 
tuting the resulting expression for \f~B into Eq. (|A143|) . 
we find 



d\n ( 



„n-n 



M 



3 k 
~2v 



y/Ko(4>)Qo\l + W 



(n-1) 



,,n— 1 



dv 



,.n—l 71 /r—1 



M dln< 



M~ 



dM 

din <p 
(A154) 



(A146) Rewriting the last equation as 



1 dv r - 



din ill 



w m = -- 



n dln< 



din cj) 
(A155) 

we note that the left-hand side tends to a positive limit 
(or to a positive infinity) due to Eq. (|A152[) , while the 
term dv n /dhx(f) tends to zero at <p — > 00 and can be 
neglected. Therefore, for large 4> we obtain 



3 k 



x/K (<j))Qo\l + w m \M 2 



dM 



dln« 



(A156) 



This relationship is sufficient for our purposes; we will 
now show that v(<j)) cannot tend to zero at <p — > 00. If 
we assume that v((f>) — ► 0, we must have 



lim -V^E = 0. 



d In <p 



M- 



(A157) 
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Using Eq. (|A151|) , we transform this condition into 



lim 

d> — >oo 



2 ( mif o),ln0 



dln< 



d\nK 
d\nd) 



(A158) 

It is now straightforward to show that the condi- 
tion (|A158p cannot be satisfied by a function Kq{4>) that 
tends to zero at cp — > oo. Since (lnifo) < for all 0, 
the function (mi^o) i n< k tends to a nonpositive constant 
or to a negative infinity at <f> — > oo. Hence, we obtain the 
bounds 



1 < 



i(lniTo) 



,ln<t> 



2-QnK ) 



< 0, < 



M4> 



2- 



d\nKo 
din 6 



< 



1 



2 

(A159) 

The derivative of a bounded function cannot have an infi- 
nite limit. Therefore the limit (|A158p cannot be infinite. 
Since the condition (| A158|) cannot be satisfied, solutions 
with v(c/>) — ► and B(cf>) — > do not exist under the 
present assumptions. 

Having shown that Ro = 1 is impossible, we assume 
Rq < 1 in the rest of this section. Let us now consider the 
admissible behavior of v(4>) at large cj). It is convenient 
to change the independent variable from <f> to x defined 
by Eq. (|A144p and to rewrite Eq. (|A14ip as 



dv 
d\ 



1 



1 



Qo 



— 3k 



Qo 



nQiw™" 1 \VKj a ""V 1 -R 

(A160) 

For an asymptotically stable solution, we need v{4>) — > 
while > 0. Therefore, dv/d<fi (and therefore also 

dv/dx) must remain negative at large 0. Let us ex- 
amine the condition under which the right-hand side of 
Eq. (|A160|) might be negative at large <f>. 

We notice that the first term in the right-hand side 
of Eq. (|A160[) contains a negative power of v multiplied 
by a nonnegative function d(K~ 1 / 2 )/d(j) and a positive 
constant. This term will diverge to positive infinity as 
v — ► unless d(K~ 1 / 2 )/d(j) tends to zero at large <f>. On 
the other hand, the second term, 



3k 



(A161) 



tends to a negative constant at large 0. Thus, dv/dx may 
become negative at large <f> only when d(K~ x l 2 }j d<\) tends 
to zero at large 4>. If K{<j>) is such that K'K~ 3 / 2 — > 0, 

then xO) = J ^/K{^)d4> diver ges at — > oo; this was 
already shown in the previous section after Eq. ()A129|) . 
Let us therefore continue the analysis under the assump- 
tions (|A129|) and Rq < 1, taking into account that 
X — > oo together with <j> — ► oo. 
Rewriting Eq. (|A142|) as 



dR 

dx 



3 k 



RVI^R^/Qo{w m + I)[l + O{v)\, (A162) 



and noting that w rn + 1 =/= 0, we immediately see that 
dR/dx — > can be realized only if w m + 1 > and either 



Ro = or Rq = 1, where i?o = hni0_»oo R{4>)- Since we 
are assuming i?o < L the only admissible value is Rq = 
0. Therefore, we now look for solutions of Eqs. (|A160|) - 
(|A162j) such that v(x) -> and R(x) -> as x -> oo (at 
the same time as — * oo). 

Computing the limit of Eq. (|A160p as — > oo and not- 
ing that dv/dx -* on asymptotically stable solutions, 
we obtain the condition 



lim 



Qo 



i 



nQi v T 



3k y Qq. 



(A163) 



The right-hand side above is a nonzero constant. There- 
fore it suffices to look for solutions v(<p) of the form 



'») = 



(A164) 



where A(0) is a new unknown function replacing v((j)). 
Solutions v(<j>) — > will be asymptotically stable if the 
general solution for A(c/f)) tends to zero as — * oo. For 
brevity, we rewrite the ansatz (|A164|1 . with the indepen- 
dent variable <fi expressed through x, as 

v(x) = [(l + A(x))W( X )}^ , (A165) 
where W(x) is a fixed function defined through 



w(x)L = 



x=x(<t>) 



3nnQi 



(A166) 



By assumption, we have IU(x) — > as x — > oo. Sub- 
stituting the ansatz (|A165|) into Eqs. (|A160|) - (|A162|1 . we 
obtain, to the leading order in A and R, 



dA 

dx 

dR 

dx 



3kVQo" {\R + A) + (n - 1) (V^r) 



= -Sk^/Qq (w m + 1) RW- — . 



(A167) 
(A168) 



Since the equation for R does not contain A, the stability 
analysis can be performed first for R(x) and then for 
A(x) assuming that R(x) ~> 0. It is convenient to replace 
the independent variable x temporarily by 



ip(x) 



W~ 



~(x) d x- 



(A169) 



Since VF(x) — + as x ~~ * 00 j the new variable ip grows 
to infinity together with x- The new equations for A{tp) 
and R{ip) are 



i 



dR _ 



?0 (u>m + 1) 



>(^) V 

(A170) 
(A171) 
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It is clear that the general solution for R(ip) tends to zero 
if w m > — 1. The general solution for A(ip) is a sum of the 
general homogeneous solution (which tends to zero) and 
an inhomogeneous solution. The inhomogeneous terms 
are proportional to R and (lU 1 ^"" 1 )) , both of which 
tend to zero at x — > oo (ip — > oo). Therefore the general 
solution for v((j)) and R((j)) is asymptotically stable under 
the current assumptions. 11 



8. Slow motion (vo = 0), marginal cases (Q(0) = 0) 

Let us now turn to the case Q(0) — 0. In this case, we 
may expand the relevant quantities near v = as follows, 



Q{v) = Qiv n [1 + 0{v)} , 
i^v) = {n-l)Q 1 v n [l + 0(v)], 
1 



C 2 s(v) = 



n- 1 
1 



n 



1 



[l + 0(v)} 
[l + 0(v)} 



(A172) 
(A173) 

(A174) 
(A175) 



where n > 2 and Q\ > 0. Using these expansions, we 
rewrite the equations of motion ([35 |) -(|36 | . in the leading 
order in u, as 



dv 
~d4> 
dR 



v K' 



V~r. 



1 V 1 - R 



-V 2 



= SkRVI - RV(n - 1) QiK(4>)- 



(A176) 

w<j,(v) 



v 



l-n/2 



(A177) 



The possible asymptotic values of equation of state pa- 
rameter w^(0) are l/(n — 1) for n > 2; in particular, 
we can have w^O) = -|, mimicking radiation, if n = 4. 
When w m = l/(n— 1), we may need to expand the term 
w m — u^f) to a higher nonvanishing order in v. For 
instance, assuming an expansion 



Q(v) = Qiv n + Q 2 v n+P [1 + 0(v)} , 
where n > 2 and p > 1 , we find 



(A178) 



Q(v) 



vQ'{v) - Q 



1 + Q(v) 
n- 1 



(n-l)Qi 

(A179) 

Let us begin by considering the possible asymptotic 
value Rq = 1 of R(<fi) at (j> — > oo; values i?o < 1 will be 
considered subsequently. In case Rq = 1, we write the 
ansatz i?(0) = 1 — -B(^) and transform Eq. (|A177p into 



dy/B 



w<t,(v)) . 
(A180) 



11 This is case 10 in Sec. lVXl 



The right-han d side of the equation above must be nega- 
tive to allow ^B(cf)) — > at <p —> oo. This cannot happen 
if w m — w,j,(0) > 0. Thus, the only possibility for the exis- 
tence of stable solutions is w^,(v) > w m for v > (which 
does not exclude w,p(0) = w m ). Under the assumption 
w$(0) > w m , Eqs. (|A176|) and (|A180|) can be rewritten 
(again keeping only the leading-order terms) as 



In [K x ' r - 



SkVQI K(<f>) 
'^TV B(<f>) 



v 2 



(A181) 



d(f> — 1 



(1 - (n - l)w m ) 



B(4>) 



v 2 



(A182) 



In case w^,(0) = w m , we assume the expansion (|A178|) 
and use Eq. (|A179|) ; then Eq. (|A182|) is replaced by 



dlnB 



(n-1) 



3/2 



B(4>) 



V 2 



-1+P 



(A183) 



As in the case u^(0) ^ w m , stable solutions are possible 
only if the right-hand side of Eq. (|A183p is negative, i.e. if 
Q 2 <0. 

We now need to analyze the solutions of the sys- 
tems (lAT8T1) - (jAT82)) and (|AT8T|) . (fATM)) by looking for 
such K{(f)) that the general solutions v(<p) and B{<f>) al- 
ways tend to zero in the two cases. 

The general solution of Eqs. (|A181^ - i|A182|) can be 
found by first noticing that 



In K^ n v) - 



IylB 



1 — (n — 1) w-n 



= 0. 



Hence we may express 



B{<t>) = Co 



l — (n — l)w r) 



(A184) 



(A185) 



where Co > is an integration constant. Then we substi- 
tute this B{4>) into Eq. (|A181[) and obtain the following 
equation for the auxiliary function u = K^-^v, 



du 



(A186) 



where we defined the auxiliary constant s and function 

F{(j)) as 



(1 



F{4>) = 3k, 



(n - 1) Co 



K^ n (<f>) > 0. 



(A187) 
(A188) 



Since by assumption w m < , the possible values of s 
are s < ^. We are now looking for functions K(<j>) such 
that both v = u/F and B oc u n ~ 2s always tend to zero 
as — > oo; in other words, we require 



<<j>) 



lim 



0. 



lim u{(j)) 







(A189) 
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for the general solution u(cj>). The general solution for 
u(4>) can be written as 



u(<f>) 



i/(i-«) 



S = 1, 



+ 

exp^- f F(0)#) , 

(A190) 

where Ci is a constant of integration. If s < 1, the power 
1/(1 — s) is positive and so the general solution u(<j>) does 
not tend to zero. If s > 1, the general solution 

tends to zero in case F(<j))d4> diverges as <f> — > oo, and 

does not tend to zero if ^F{<f>)d<f> converges. Therefore, 
the only possibility for a stable solution is s > 1 and 
J F{(j))d(j) — > oo as <fi — > oo, or equivalently 

w m >-^—^\ hm / K l l n (<f))d<j> = oo. (A191) 



1 



un- 



it remains to examine the condition u(<j>)/ F((f>) 
der these assumptions. 

Since we already have u — > 0, the condition u/F — > 
holds if -F(0) approaches a nonzero constant or infinity 
as </> — ► oo. However, if 



lim F(<j>) = lim K 1/n (4>) = 0, 



(A192) 



the condition it/i* 1 — > is a nontrivial additional con- 
straint on the function K ((/)). This constraint can be 
expressed as a condition on K(<f>) as follows. We find 
from Eq. (|A190|) that 



u 

— oc 
F 



F^ 1 J* F((j>)d(j> , s>l, 

exp (-lnF- ^ F{(j))d(^j , s = 1. 

The condition — > is then equivalent to 



(A193) 



lim F s 

(f> — >oo 

lim In F 4 



oo, s > 1; (A194) 



oo, s = 1. (A195) 



We note that the left-hand sides in Eqs. (|A194[> — (|A195[> 
depend monotonically on the growth of F((f>); more pre- 
cisely, the terms under the limits become larger when 
we choose a faster-growing or slower-decaying function 
F(4>). Thus, it is clear that the conditions (|A194j) - (|A195|) 
will hold if F((f>) decays sufficiently slowly as <j> — * oo 
(or grows, but this case was already considered). With 
some choices of F(<fr) = F (</>), the limits in Eqs. (|A194|I 
(IA195P will be finite nonzero constants. We can easily 
determine such i*b (</>), 

F {<f>) oc (j)~ 1/s , s > 1, (j>-> oo. (A196) 

Hence, the limits (|A194|) — (|A195|> will be infinite when 
F{4>) decays slower than The corresponding con- 

dition for K((f)) can be written as 



lim (t> n/s K(4>) = 



(A197) 



We can make this argument more rigorous by assuming 
the ansatz 



(A198) 



where K\(4>) > is an auxiliary function. Note that the 
function F{<f>) is related to K{<f>) by Eq. (|A188|) . which 
contains an arbitrary integration constant Co > 0. Thus 
we may write 

Ftf) = C 1( f>-^ s Kl /n (4>), (A199) 
where C\ > is an arbitrary constant. If 



lim Ki{4>) = oo, 

<fi — >oo 



(A200) 



it means that K\{<p) is larger than any constant at suffi- 
ciently large cj>. Then we obtain lower bounds (for arbi- 
trary constant C*2 > 0) 



F(4>)d4> > C 2 ^ 1/s , s > 1; (A201) 
r<t> 

/ F{<j))d<j) > C 2 ln0, a = l, (A202) 



and the conditions (|A194jK(|A195[) hold. On the other 
hand, if 



lim Kxty) = K[ 0) < oo, 

(p—>-OG 



we find 



1-1/ s 



s > 1; 



(A203) 



(A204) 



r<t> 

j F(<p)d<p& C 2 ln<p, s = l, (A205) 

where C2 > is an arbitrary constant. In that case, 
the conditions (|A194[) — (|A195[) cannot hold for arbitrary 
C<z- Therefore, the condition (j A197|) is necessary and 
sufficient for Eqs. (|A194|) - (fAT95l) to hold. 

We conclude that an asymptotically stable solution ex- 
ists for vo — 0, Q(0) = with the expansion (|A1 781) . 
when Rq = 1, w m < ^ry, and the conditions (|A191[) 
and (|A197|) hold. We note that for n = 2 the condition 
w m < Tjzj contradicts the first condition in Eq. (|Al91jl . 
so admissible solutions exist only for n > 2. 12 

The remaining case requires the analysis of 
Eqs. (|A18ip — (|A183|> . The general solution of these 
equations cannot be obtained in closed form; however, 
we only need to analyze the asymptotic behavior at 
(j) — > 00. So we will estimate the relative magnitude 



12 This is case 6 in Sec. W~K\ 



23 



of different terms in these equations. Let us rewrite 
Eqs. (|M8T|) - (rM83)) as 



dlnv 1 K' ~ K 1 / 2 r 



■i/2-l 



-Q 2 K l l 2 v n ' 2 ~ l+ P, 



where the auxiliary positive constants 



Qi 



i 



(n-l) 



3/2 



(A206) 
(A207) 

(A208) 



were introduced for brevity. (Positivity of these constants 
is clearly necessary for the existence of asymptotically 
stable solutions.) Suppose that v((f>) and B{(f)) are de- 
caying solutions of Eqs. (|A206|) - (IA207p . and let us com- 
pare the magnitude of the terms in the right-hand side of 
Eq. (|A206[) in the limit <fi — ► oo. There are only three pos- 
sibilities: the first term dominates; the two terms have 
the same order; or the second term dominates. In other 
words, the limit of the ratio of the second term to the 
first, 



Q lK 3/2 v n/2-l 

q = hm == . 

4,^00 K'VB 



(A209) 



must be either zero, or finite but nonzero, or infinite. The 
value of q must be the same for every decaying solution 
{v((f>), B(ip)} except perhaps for a discrete subset of solu- 
tions, which we may ignore for the purposes of stability 
analysis. In each of the three cases, Eqs. (|A206I) — (|A207[) 
are simplified and become amenable to asymptotic anal- 
ysis in the limit (f> — > oo. We will now consider these 
three possible values of q in turn. 
If q = 0, we have at large cf> 



1 K' ~ K 1 / 2 



n K Jb 



n/2-l 



(A210) 



and thus only the first term is left in Eq. (|A206|1 . 

dlnv IK' n 

— — « =4> v(xK~ 1/n . A211) 

dqy n K 

Decaying solutions have v((f) — > 0; so a necessary condi- 
tion is K- 1 / n {4>) -> at <j) -> oo. With this v(4>), the 
condition (|A210|) becomes 



(W«Y» % 



(A212) 



However, this condition cannot be satisfied, since the left- 
hand side tends to zero at large <p while the right-hand 
side tends to infinity because B — > 0. Thus, decaying 
solutions v((j)), B(cj)) are impossible with q = 0. 

If q ^ and |<7| < oo, we consider q as an unknown 
constant that possibly depends on the solutions v((f>) and 



B{4>). At large (j>, we have 



K 



dlnv K' f I 



Therefore, 



v{(j>) oc K- q - 1/n ((j)) -» 



(A213) 
(A214) 

(A215) 



since we need a decaying solution v((j>). With this v (<fi), 
Eq. (|A213|) yields 



r- Qi ifi-(«/2-i)9+lA 

V ^ w — tF, 

q K 



(A216) 



For a decaying solution B((f>) — ► 0, we thus must have 



1 d 

^ "77 

B dd> 



K (n/2-l)q-l/r, 
and in particular 



oo as (f> — > oo, (A217) 



n \ 1 . 
? -i)g--^0. 

2 / n 



(A218) 



Substituting the expressions for v{4>) and y/Bl^) into 
Eq. (TA"2"07l) . w e nnd 



dy/B 



d(j) q 



K' 







1 - 


(§-') 



1 ifiT" 

n 



if' 2 



(A219) 



This is now a closed equation for K (</>), which we may 
rewrite as 



if y 
f> 1 



if if" 
if' 2 



[(§-0 








n 



(A220) 



Due to the conditions ljA215|) . (|A218j) . and since p > 0, 
the right-hand side above tends to a nonzero limit as 
<f> — > oo, namely 



if V 
if 



(?-') 



(A221) 



It follows that the only admissible form of the function 
if (0) is 



if(»oc0 Q , 0^oo. 



(A222) 
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However, this expression does not satisfy Eq. (|A217[) . 
Therefore, asymptotically stable solutions are impossi- 
ble. 

In the last case, q — oo, we may disregard the first 
term in Eq. (|A206[) and obtain 



dlnv ~ K 1 ' 2 



,«/2-l 



(A223) 



Then we can rewrite Eq. (|A207[) as 



Q 2 d 

pQi d( t> 



(A224) 



This relationship between B and v can be integrated and 
yields 



Vb~ = exp 



pQi 



(A225) 



where C\ is a constant of integration. It follows that it 
is impossible to find simultaneously decaying solutions 
v(4>) — ► and B{4>) — * at <j> — ► oo. 

This concludes the consideration of the case Rq = 1 
and w m — -p^j , in which case there are no asymptotically 
stable solutions. 

We now turn to the analysis of the case Rq < 1- We 
first note that the leading terms of Eq. (|A176[) do not con- 
tain R when R — > Rq < 1. Therefore the stability analy- 
sis can be performed for R{4>) and v(<p) separately. Using 
the ansatz R{<j>) = Ro + B(</>) and assuming a fixed solu- 
tion v ((/)), we find that the right-hand side of Eq. (|A177|) 
is independent of B{4>) if < Ro < 1. Therefore, general 
solutions B((f>) will not approach zero in case Rq ^ 0. It 
remains to look for asymptotically stable solutions v(<p) 
and R{4>) in case Rq = 0. 

In case Rq = 0, we begin by analyzing the asymptotic 
behavior of v((f>). Rewriting Eq. (|A176[) as 



du 



u 2 K 



l/n 



(A226) 



we find the approximate general solutions (valid only for 
large <p) 



i(<f>) = exp 



[■<t> 

J O 



«(0) 



K l / n d<t> 



2, (A227) 



n > 2, (A228) 



where 4>q is an integration constant. The general solution 
v = K~ 1 l n u should tend to zero as — > oo. We note that 
Eq. (|A226[) is similar to Eq. (|A186|) after the replacements 

Ffe) = ^SbfV-fe), , s » (A229) 
y/n—l 2 



Therefore, we may use the conclusion obtained after 
Eq. (|A190j) . with the caveat that F(<fr) is presently re- 
lated to K(cj>) uniquely, without an arbitrary proportion- 
ality factor. This was used to exclude the boundary 
case (|A203|) . which is presently still allowed. Thus the 
condition (|A197|) obtained above, 



lira (/> n/s K((f)) = lim cj> 2 K{cj)) = oo, 



(A230) 



is now merely a sufficient condition for the stability of 
the general solution v (</>). In the boundary case, 



lim (j> 2 K(4>) = Kq, < Kq < oo, 



(A231) 



we find 



v{4>) oc exp ^1 - 3k VQi^oJ ln^ , n = 2, (A232) 
v((f>) w const, n > 2. (A233) 

Thus, the case (|A23ip yields a stable solution for v((f>) 
when n = 2 and 3k\/Q±Kq > 1. (The possibility 
Sk^/QiKq — 1 is unphysical because it requires an in- 
finitely precise fine-tuning of the parameters in the field 
Lagrangian.) Thus a sharp condition for the asymptotic 
stability of v((f>) is 



K{cj>)> 



1 



at 



oo, n = 2; 



(A234) 



lim 4> 2 K{(j)) = oo, n > 2. (A235) 



A weaker necessary condition is 

K l / n {<j>)d(j) = 



S 



(A236) 



Let us now consider the stability of the general solution 
for R(4>). It follows from Eq. (jA177|) that 



d\nB 



This equation integrates to 



(A237) 



5(0) = Bq exp 



-const ■ 



(A238) 

where Bq is an integration constant. The general solu- 
tion for B(4>) will tend to zero as long as the integral in 
Eq. (|A238j) diverges to a positive infinity at oo, 



{w m - Wj,(v)) K^v 2 l d(j) 



(A239) 



A necessary condition for that is w m > — ^r. Precise 
constraints on K(<f>) for Eq. (|A239|) can be obtained by 
considering the cases n = 2, n ^ 2, w m = -^zj, and 
w m > separately. 
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If w m > —^j , the condition (|A239[) holds when 



K^v^~ l d(j) = oo. (A240) 



If n = 2, the above integral diverges due to the necessary 
condition (|A236|) . If n > 2, we use the solution (|A228|) . 
where u = if 1 /™?;, to obtain 



= / # 



Temporarily introducing the auxiliary function 



m = / k^mm, 



K 1/n ((f)). 
(A241) 

(A242) 



we note that lim^oo I ((f) = oo by Eq. (|A236|) . Therefore 
we express Eq. (|A241|) through I ((f) and obt am 



(A243) 



lim \nl((j)) + const = oo. 

(j) — * oo 



Therefore, the general solution B(4>) tends to zero with 
w m > — zr f° r an Y n>2 under the condition (|A236[) . 

When w m = ^zj, it follows from Eq. (|A178|) that the 
integrand in Eq. (|A238[) acquires an additional factor pro- 
portional to v p , where p > 1. Therefore the general solu- 
tion for B((f) will tend to zero only if 



K l/2 v n/2-l+p d(j) 



OO, 



(A244) 



where we need to substitute v((f) = K 1 / n u(^) and u((f) 
as given by Eqs. (|A227p - (|A228|l . 

Consider first the case n — 2; we will now show that 
the condition (|A244[) is incompatible with the earlier con- 
dition (|A234[) . Using the solution (|A227|) , we can rewrite 
the condition (|A244|) as 



,#(i-j>)/2 



exp 



r<t> 

J <f> 



(A245) 



By the condition (|A234I) . we have 

jf(i-p)/a < const . ^-^ ( A 246) 

/ if 1/2 (0i)d0i > V^o In + const. (A247) 

J 4> 



Therefore the integral in Eq. (|A245p is bounded from 
above by 



const 



const 



" < oo, 
(A248) 



where we temporarily denoted 

a = (^Ky/QiKo - l) p > 0, (A249) 
and so the condition (|A245p cannot hold. 



It remains to consider the case w„ 



and n > 2. 



Using Eq. (|A228j) . we rewrite the condition (|A244j) as 



K^((f) 



<t>0 



According to Eq. (|A235j) . we must have 

K(t/>) > C^ 2 



= oo. 
(A250) 

(A251) 



for any Co > at large enough d>; thus K((f) should decay 
slower than (j>~ 2 . However, it is straightforward to verify 
that a power-law behavior 



K((f) oc (f 



-2+S 



oo, S > 0, 



(A252) 



yields a convergent integral in Eq. (|A250|) . Therefore, 
the only possibility of having an asymptotically stable 
solution is to choose K(<j>) such that it decays slower than 
<p~ 2 but faster than (f)~ 2+s for any 5 > 0. An example of 
an admissible choice of K(<f) is 



K{4>) oc d>- 2 (lnc 



a > 0. 



(A253) 



With this K ((p), we obtain the following asymptotic es- 
timate at large <fi, 



iv 1/,l (0i)#i oc const • (j) 1 - 2 ' 11 (ln</>) Q/n , (A254) 



<t>0 



and so the integral (|A250[) becomes, after some algebra, 

P 



const 



(j)- 1 (lnc 



i-2 rlrh — 



if, 



71-2 



< 1. 



(A255) 

Since the convergence of the integral in Eq. (|A250[) mono- 
tonically depends on the growth properties of the func- 
tion K ((/)), it is clear that the condition (|A250[) will also 
hold for functions K(cf) satisfying Eq. (|A235|) but grow- 
ing slower than those given in Eq. (IA253|) . However, the 
condition (|A250[) may not hold for K ((f) growing faster 
than those in Eq. (|A253|) . 

To investigate the admissible class of functions K ((f) 
more precisely, let us use the ansatz 



K((f) = (f>- 2 K (4>), 



(A256) 



where Ko(<f) is a function growing slower than any power 
of 4>. Then we have an asymptotic estimate (for n > 2) 



K^^dfa » const • (f- 2ln (K Q ((f)) lln , (A257) 
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and we can rewrite Eq. (|A244p as 



const ■ 



(f)- 1 [K ((j>)]~*=* d(f> = oo. (A258) 



Substituting K = (j) 2 K into Eq. (|A258p . we find that the 
conditions (| A235|) and (|A250p are equivalent to 

/oo 
[X d(/> = OO. 

(A259) 

The condition ()A235|1 guarantees the stability of v(<f>), 
while Eq. (|A250[) guarantees the stability of B(<p). There- 
fore, Eq. (|A259[) is a sharp (necessary and sufficient) con- 
dition for the stability of the solution {v, B}. 

A sufficient (but not a necessary) condition for the di- 
vergence of the integral in Eq. (|A258|> is 

lim (ln0)~^ K (4>) < oo. (A260) 
The corresponding sufficient condition for K (<fi) is 



lim cj) 2 K((j)) = oo, lim (f) 2 (m<£) * K{4>) < oo. 

(A261) 

13 This is case 6 in Sec. IVAl 

14 This is case 7 in Sec. IVAl 

15 This is case 11 in Sec. lVXl 

16 This is case 12 in Sec. N~K\ 



The sharp condition (|A259[) cannot be restated in terms 
of the asymptotic behavior of K(<f>) at 4> — > oo, but of 
course one can check whether Eq. (|A259p holds for a given 
K((fi). The condition (j A259|) specifies a rather narrow 
class of functions; however, we strive for generality and 
avoid prejudice regarding the possible Lagrangians. 



In this section we have shown that asymptotically sta- 
ble solutions exist with vo = and Q(Q) = only in the 
following cases: (a) Asymptotic value Rq = 1. Expan- 
sion (|A172[) holds with Q\ > 0, determining the value of 
n, which should be n > 2; — 2^3 < w m < accord- 
ing to Eq. (|A191|) ; a nd K[ 4>) satisfies Eq. (|A197p . where 
s is defined by Eq. (|A187p . 13 There are no stable solu- 
tions when w m — —^j and expansion (|A178p holds, (b) 
Asymptotic value Rq — 0. Expansion (|A172p holds with 
Qi > 0, determining the value of n > 2; either n = 2, 
w m >l, and K ((f)) satisfies Eq. (|A234|) 14 or E q (jM35]l . 15 
or n > 2, w m = ^j, and K{<j>) satisfies Eq. (|A259p . 16 
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